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ABSTRACT 

Any  slowly  varying  linear  space- variant  system  can,  in 
principle,  be  represented  holographically  by  spatially  samp¬ 
ling  the  input  plane  and  multiplexing  the  respective  system 
transfer  functions.  A  scheme  reported  earlier  for  implement¬ 
ing  this  technique  makes  use  of  phase  diffusers  in  the  ref¬ 
erence  beam  paths  to  encode  sequentially  recorded  holograms. 
However, to  minimize  the  cross  talk  between  the  holograms  upon 
playback  the  diffusers  should  have  good  correlation  properties. 
In  this  report  extensive  computer  simulations  to  evaluate  the 
correlation  properties  of  a  family  of  binary  phase  codes  are 
conducted.  An  alternative  multiplexing  technique  in  which 
the  transfer  functions  are  sampled  in  the  Fourier  plane  to  gen¬ 
erate  a  composite  hologram  is  also  described.  In  this  tech¬ 
nique  the  samples  of  the  transfer  functions  are  placed  in 
nonoverlapping  regions  and  hence  there  will  be  no  crosstalk 
upon  playback.  However  multiple  copies  of  the  input  function 
are  required  during  the  playback  step.  The  results  of  prelim¬ 
inary  experiments  conducted  to  evaluate  this  approach  for 
space-variant  system  representation  are  presented  including 
the  verification  of  coherent  addition  using  computer  multi¬ 
plexed  holograms. 
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CHAPTER  1 


INTRODUCTION 

A  method  for  holographically  representing  any  band- 
limited  space-variant  system  using  a  sampling  technique 
has  been  described  [1,2].  This  method  requires  the  se¬ 
quential  multiplexing  of  a  number  of  holograms  of  the  system 
transfer  functions  in  a  single  recording  medium.  As  a 
result  of  this  multiplexing,  many  crosstalk  terms  are 
generated  upon  playback  in  addition  to  the  required  system 
response  terms.  Several  crosstalk  suppressing  techniques, 
using  properties  such  as  extinction  angle  effects  of  volume 
holograms  [3] ,  and  the  correlation  properties  of  phase  codes 
when  used  in  the  reference  beam  paths  [4,5]  have  been  suggested 
for  implementing  this  scheme.  Experimental  results  using 
ground-glass  diffusers  and  binary  amplitude-coded  diffusers 
have  also  been  reported  [6] .  Analytical  studies  to  model 
the  characteristics  of  various  diffusers  have  been  carried 
out  [7] .  Use  of  randomly  generated  binary  amplitude  diffusers 
with  computer  multiplexed  holograms  has  also  been  studied  [8] . 
Preliminary  studies  on  diffusers  based  on  the  known  correla¬ 
tion  properties  of  the  Gold  codes  used  in  spread  spectrum 
communication  systems  have  been  carried  out  [9] . 

The  work  presented  in  this  report  consists  mainly  of 
two  parts.  In  Chapter  2  results  of  extensive  computer 
simulations  to  study  the  auto  and  cross-correlation  properties 
of  Gold  codes  of  various  lengths  under  different  conditions. 


2 

along  with  the  computer  simulated  outputs  when  these  codes 


are  used  as  phase  diffusers  in  multiplex  holography  are 
presented.  In  Chapters  3  and  4  an  alternate  method  of  multi¬ 
plexing  the  transfer  functions  using  a  sampling  technique, 
along  with  preliminary  experimental  results  using  computer 
multiplexed  holograms,  is  presented. 

1.1.  Sampling  Theorem  for  Space-Variant  Systems  [1] 

The  output  g(x)  of  a  linear  system  due  to  an  input  f(£)  is 
given  by  the  superposition  integral 

g(x)  =  S  [f  (?)  ] 

00 

=  /  f(e)h(x-£r€)d£  (1-1) 

—00 

where  S[*]  is  the  linear  system  operator.  The  system  line- 
spread  function  h(x-£»£)  is  the  system  response  to  an  input 
Dirac  delta,  [10] 


h(x-5,5>  =  S  [5  (x-?)  ]  (1-2) 

Now,  Fourier  transforming  the  Eqn.  (1-1)  we  obtain 

G(fx)  =  Fx(g(x)) 

00 

=  /  f  (£)Fx[h(x,C)  ]exp(-j2*fxC)d£ 

—00 

=  F  r  F  [f  ( 5)  h  (x  ,  5)  ]  I  _  (1-3) 

s  x  I  v=  f 

X 

Where  v  and  fx  are  the  frequency  variables  associated 
with  C  and  x  respectively. 


Defining  the  system's  spatial  transfer  function  as 


W55  *  Fxth(x,C)]/  (1-4) 

Equation  (1-3)  may  be  rewritten  as 

G<fx)  *  F^t  f  (OH^  (^  ,  C)  1  d-5) 

V"fx 

If  f  (£)  and  h(x,£)  are  band-limited  in  v  and  have  respective 
band  widths  of  2w^  and  2w^,  then  the  total  band  width  of 
their  product  is  given  by 

2w  *  2wf  +  2wy.  (1-6) 


Equivalently, 

g(x)  =  I  f  (£  )h(x-E  ,£  )  *  Sine  (2wx ) ,  (1-9) 

_  n  n  n 

n 

where 

Sine  (fx)  =  F [Rect  (x) ] . 


Thus  when  the  input  function  f(£)  and  the  line  spread  func¬ 
tion  h(x,£)  are  band  limited,  the  output  g(x)  of  the  system 
can  be  computed  exactly  by  sampling  the  product  of  the 


functions  f ( C )  and  h(x,5)  at  intervals  of  l/2w  and  passing 
the  sum  of  these  sampled  products  through  a  suitable  low  pass 
filter. 

Although  Eqn.  (1-9)  implies  a  countably  infinite  number 
of  samples  of  the  product  of  the  spatially-varying  system 
response  h(x-£,£)  and  the  input  function  f(£),  in  practice, 
if  fi(£)  is  essentially  zero  outside  the  interval  |£|<a  and 
if  the  spectrum  of  f(£)h(x,£)  is  essentially  zero  outside 
of  the  interval  |v|  <_  w,  then  the  required  number  of  samples 
for  a  good  approximation  is  given  by  the  space-band  width 
product 

N  =  4wa.  (1-10) 

Two  possible  schemes  for  implementing  this  sampled  system 
representation  will  be  discussed  in  the  following  sections. 

1.2.  Representation  of  Space-Variant  Systems  Using  Phase 
Coded  Reference-Beams 

A  scheme  for  coherently  representing  a  space-variant 
system  using  the  sampling  technique  described  in  Section 
1.1  is  shown  in  Figures  (1-1)  and  (1-2),  [4].  During  the 
recording  step  the  space  variant  system  is  sequentially 
sampled  in  the  input  plane  at  N  points  denoted  by  i^  through 
iN  to  generate  the  spread  functions  h^  through  h^.  The 
corresponding  reference  beam  diffuser  functions  are  de¬ 
noted  as  r^  through  r,j.  After  Fourier  transformation  by 
lenses  and  the  amplitude  transmittance  t  of  the 


Playback  scheme  for  the  holographic  representation 
of  space-variant  system. 
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hologram  is  given  by 

N  2 

t «  i  i  (h.  +  r.) r  ,  (1-id 

i*l  1  1 

where  h.  and  R.  are  the  Fourier  transforms  of  the  functions 

i  i 

hi  and  r^  respectively. 

During  playback  the  input  function  f(£)  is  spatially 

sampled  at  the  sample  locations  of  the  original  reference 

array  to  produce  the  sampled  inputs  through  where 

til 

is  the  sampled  value  of  the  input  function  at  the  i 
location.  Then  the  reconstructed  wavefront  to  the  right 
of  the  hologram  is  given  by 


The  term  R. 

i 


H.  +  R. 

i  i 


may  be  expanded  as 


Ri|Hi  +  Ri|2  =  RiR*Hi  +  +  RiRiHi  +  RiRiRi 


(1-13) 


where  *  represents  the  complex  conjugate  operator. 

* 

Out  of  these  components  only  the  term  Rj_R^H^  Is  diffracted 
by  the  hologram  in  the  direction  of  the  output  plane  as  a 
result  of  the  offset  in  the  reference  beam  path.  Hence 
the  output  after  Fourier  transformation  by  lens  is 


gggHgg  *  -  - -  - 
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given  by 


g'  (x)  - 


N  N  N 

Z  sihi  *  +  I  l  sihi  *  <ri*r  J' 

i-1  1  1  1  1  i=l  j-1  3  3  1  3 

^  (1-14) 


where  ★  denotes  Correlation  and  *  denotes  Convolution. 
The  desired  output  g(x)  is 

N 

g ( x)  =  l  s.h.  •  (1-15) 

i-1  1  1 


Thus  for  perfect  reconstruction  of  the  output  g(x)  the  dif¬ 
fusers  in  the  reference  beam  path  should  have  the  follow¬ 
ing  characteristics. 


ri(x)  *  r  (x)  -  5(x)  for  all  i, 
r\(x)  k  r  j  (x)  *  0  for  all  i  and  j,  i^j 


(1-16) 


This  implies  that  we  need  a  set  of  codes  for  the  fabrication 
of  diffusers  having  delta  like  autocorrelations  and  zero 
cross-correlations.  However  such  a  set  of  ideal  codes  does 
not  exist.  In  Chapter  2  the  properties  of  a  family  of  dif¬ 
fuser  function  derived  from  a  set  of  codes  known  as  Gold 
Codes  and  previously  used  in  spread  spectrum  communication 
systems  are  studied  through  a  series  of  computer  simulations. 


1.3  Representation  of  Space-Variant  Systems  Using  a  Sampled 
Input/Sampled  Transfer  Function  Approach 
The  sampling  theorem  of  Section  1.1  requires  the  multi¬ 
plexing  of  a  number  of  system  transfer  functions  on  a  single 


t 
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hologram/  with  the  capability  of  independent  access  to  each 

*  of  the  transfer  functions.  In  Section  1.2  a  method  for 
achieving  this  independent  access  requirement  by  using  dif¬ 
fusers  in  the  reference  beam  paths  was  discussed.  However 

*  when  the  system  line  spread  function  h(x/5)  is  space  limited 
in  X/  an  alternate  approach  in  which  the  transfer  functions 

are  sampled  in  the  frequency  plane  to  multiplex  a  number 

I 

of  transfer  functions  on  a  single  recording  medium  may  be 
employed.  This  is  in  addition  to  the  sampling  in  the  in¬ 
put  plane  as  required  by  the  sampling  theorem  of  Section 
i 

1.1.  This  technique  is  described  in  detail  in  Chapter  3. 
Preliminary  experimental  results  using  this  method  are 
presented  in  Chapter  4. 

> 

I 

> 

t 

i 
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PHASE  CODED  REFERENCE  BEAM  APPROACH 

In  section  1.2  a  previously  developed  scheme  for 
representing  a  space-variant  system  using  phase  coded  ref¬ 
erence  beams  was  presented.  It  was  also  shown  [Eqn.  (1-16)] 
that  for  ideal  playback  of  the  system  responses  the  codes 
used  in  the  multiplexing  step  should  have  delta-like  auto¬ 
correlations  and  zero  cross-correlations.  In  this  chapter 
the  properties  of  a  set  of  codes  known  as  the  Gold  codes 
are  evaluated  for  use  as  phase  diffusers  to  multiplex  a 
number  of  transfer  functions  in  a  single  recording  medium. 

In  section  2.1  a  method  for  generating  a  set  of  Gold  codes 
is  described.  In  section  2.2  the  correlation  properties 
of  the  Gold  codes  of  various  lengths  are  evaluated  through 
a  series  of  computer  programs.  The  output  of  a  system  us¬ 
ing  the  Gold  codes  for  multiplexing  the  transfer  functions 
of  a  space-variant  system  are  also  simulated  on  the  com¬ 
puter  . 

2.1  Generation  of  Gold  Codes 

An  analytical  technique  for  constructing  a  large 
family  of  codes  having  uniformly  low  cross- correlations 
has  been  described  by  Gold  [12] .  The  following  steps  des¬ 
cribe  the  technique  for  generating  a  set  of  the  Gold  codes. 
The  method  is  illustrated  by  an  example  at  the  end  of  this 
section. 


* 
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1.  Find  the  order  of  the  primitive  polynomial  required 
by  using  the  equation 

L  -  2n  -  1  (2-1) 

where  n  is  the  order  of  the  polynomial  and  L  is  the  length 
of  the  code.  The  primitive  polynomials  for  each  degree 
have  been  tabulated  [13]  . 

2.  Select  a  pair  of  preferred  polynomials  f^(x)  and  f2(x) 
that  result  in  sequences  with  low  cross-correlations.  (This 
step  is  explained  in  detail  while  describing  the  method  with 
an  example.) 

3.  Find  the  product  of  the  polynomials  obtained  in  the 
previous  step  to  obtain 

f'(x)  -  fL(x)  f2(x)  (2-2) 

4.  Convert  the  coefficients  of  the  powers  of  x  in  the 
polynomials  f'(x)  to  modulo  2  to  obtain  f(x). 

5.  Enter  these  coefficients  as  connections  of  a  2n  stage 
shift  register.  Here  0  denotes  no  connection  and  1  denotes 
the  presence  of  a  connection. 

6.  Select  a  2n  bit  binary  seed  as  input  to  the  shift 
register.  The  output  of  the  shift  register  is  a  sequence 
of  period  L  and  represents  a  Gold  code. 

7.  To  generate  another  member  in  the  set  select  a  seed 
that  is  not  a  2n  bit  segment  of  the  codes  already  generated. 
Repeating  this  step  a  total  of  (2n  +1)  sequences  each  of 
length  (2n  - 


1)  may  be  generated. 
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It  has  been  shown  that  the  cross- correlations  cc(t) 

between  any  pair  of  these  sequences  obey  the  inequality 

n+1 


cc (t)  |  < 


for  n  odd> 


n+2 

-~T~ 


(2-3) 


2  +1  for  n  even  and  n^O  mod  4, 


where  the  cross-correlation  cc(t)  between  two  codes  is  de¬ 
fined  as 


cc(t)  *  (number  of  agreements  -  number  of 
disagreements) , 


(2-4) 


when  two  codes  with  a  displacement  of  t  between  each  other 
are  compared.  The  autocorrelation  Ac(t)  for  t  ^  0  also 
obeys  the  inequality  of  Eqn.  (2-3)  and  when  t  *  0  is  equal 
to  the  length  of  the  sequence. 


Ac ( 0 )  =  2n  -1 


(2-5) 


Carter  [14]  has  described  a  method  for  generating  a  family 
of  codes  for  n  =  0  mod  4,  i.e.,  for  n  =  4,8,12,16  etc. 


2.1.1  Example  of  Generation  of  a  Set  of  511  3it  Gold  Codes: 
method  just  described  for  generating  a  family  of  Gold  codes 
is  illustrated  here  by  an  example.  In  this  example  a  set 
of  codes,  each  with  a  length  equal  to  511  bits,  is  generated 
The  computer  program  used  for  this  example  along  with  a  set 
of  9  codes  of  511  bits  each  generated  by  the  program  are 


The 
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given  in  Appendix  A.  The  following  calculations  correspond 
to  the  steps  described  earlier  for  the  generation  of  the 
codes. 

1.  The  order  of  the  primitive  polynomial  n,  is  obtained 
from  Eqn.  (2-1)  as 

511  -  2n  -1 
n  =  9. 


A  table  of  the  primitive  polynomials  of  order  9  is  given 
in  Table  (2-1)  [13]. 

In  this  table  the  polynomials  are  listed  in  octal 
notation.  For  example, 

1021  corresponds  to  001,000,010,001  and  represents 
the  polynomial  1.  x9  +  0.x8  +  0.x7  +  0.x6  +  0.x5 
+  1.x4  +  0.x3  +  0.x2  +  0.x1  +  1  i.e.,  1  +  x4  +  x9 . 
The  interpretation  of  the  numbers  in  the  first  column  is 
as  follows. 


Let  a  be  the  root  of  the  polynomial  1021.  Then 
the  number  17  in  the  first  column  of  polynomial  1333  repre¬ 
sents  that  a17,a17 ' 2l,a17’ 2  ,a17*2  ,a17*2  ,a^7*2  , 

0(17.26-5U)jC,(17.27-511)i  and  a < 17 . 28-Sll)  are  the  roots 

of  the  polynomial  1333.  Here  the  powers  of  a  are  taken  as 


modulo  511. 


1 

1021 

23 

1751 

53 

1225 

3 

1131 

25 

1743 

55 

1275 

5 

1461 

27 

1617 

73 

0013 

7 

1231 

29 

1553 

75 

1773 

9 

1423 

35 

1401 

77 

1511 

11 

1055 

37 

1157 

83 

1425 

13 

1167 

39 

1715 

85 

1267 

15 

1541 

41 

1563 

17 

1333 

43 

1713 

19 

1605 

45 

1175 

21 

1027 

51 

1725 

2.  Now  we  need  to  select  a  pair  of  polynomials  with  low 
cross  correlations.  The  first  polynomial  is  chosen  to  be 
1021.  Thus 


f^ (x)  »  1  +  x4  +  x9. 


Now  the  second  polynomial  f2(x)  must  be  chosen  such  that 
it  has  the  roots  n+^ 


(2  *  +1) 


a 

ioe.,  the  roots  of  f2(x)  must  be  a33.  Now  we  refer  to  the 
1st  column  in  the  Table  (2-1)  to  find  the  number  33.  But 
33  is  not  listed  in  the  table.  However  33  may  be  written 

C 

as  544  modulo  511  and  544  *  17. 2  .  Thus  the  required  poly¬ 
nomial  f2(x)  is  the  one  which  has  an  entry  17  in  the  first 
column.  The  corresponding  polynomial  is  1333  in  octal  repre¬ 
sentation.  Thus, 


f  2  (x) 


1  +  x  +  xJ  +  x  +  x 


7  9 

+  X  +  X  . 


3.  The  product  of  f^(x)  and  f2(x)  gives 

f  '(x)  =  fx(x)  •  f2(x) 

*  1  +  x  +  x3  +  2x4  +  x5  +  x6  +  2x7  +  x8 
+  2x8  +  2x^  +  x^  +  x^-2  +  2x^-3  +  x^-5 
+  x16  +  x18. 


4.  When  the  coefficient  of  powers  of  x  in  f^(x)  is  taken 
as  modulo  2  we  obtain 


< 


I 


I 
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f  ( ^  j  s  ^  ^  3  ^  ^  ^  ^  6  ^  ^  1 1  ^  ^  1 2  ^  ^  1 3 

+  x16  +  x18  . 

5.  The  coefficients  of  the  powers  of  x  in  the  above  poly¬ 
nomial  are  entered  as  the  shift  register  connections  to  the 
computer  program  given  in  Appendix  A.  Thus  the  shift  register 
connections  read  through  a  data  card  in  the  program  are 
101100110010110101  starting  from  the  highest  power  of  x  and 
ignoring  the  constant  1. 

6 .  The  computer  program  selects  the  seed  for  the  first  code 
as  000  000  000  000  000  001  and  generates  the  first  code 

in  the  set.  This  seed  is  then  incremented  by  1  and  the  new 
seed  is  checked  to  verify  whether  it  is  a  segment  of  the 
code  already  generated.  If  so  the  seed  is  rejected  and  a 
new  seed  is  obtained  by  incrementing  the  value  again  by  1. 

When  a  seed  which  is  not  a  segment  of  the  previously  gen¬ 
erated  code  is  found  then  the  program  computes  the  next 
member  in  the  set  of  codes.  The  program  is  written  to  gen¬ 
erate  a  maximum  of  25  sequences  out  of  the  possible  513 
sequences  that  exist  for  this  order  of  the  polynomial. 

The  auto-correlation  and  the  cross-correlations  of 
these  511  bit  codes  are  given  by  the  Eqns.  (2-3)  and  (2-5). 

Ac (0)  =  511, 

|Ac  (t)  |  <_  33  for  t  ^  0 
and  |  cc  ( t )  |  _<  33. 


17 


A  program  similar  to  the  one  just  described  for  generating 
127  bit  Gold  codes  is  given  in  reference  [15]  .  A  set  of 
nine  codes  each  of  length  127  bits  generated  by  this  program 
is  given  in  Table  (2-2)  .  In  this  table  the  binary  elements 
in  the  code  are  entered  as  0's  and  2's.  In  the  next  section 
the  correlation  properties  of  these  codes  are  evaluated 
through  a  number  of  computer  programs. 

2.2  Evaluation  of  Gold  Codes  as  Phase  Diffusers  in  Multiplex 
Holography 

In  this  section  the  results  of  evaluation  of  auto  and 
cross-correlation  properties  of  Gold  codes  of  different  lengths 
are  presented.  Also  the  computer  simulated  outputs  of  a  space- 
variant  processor  implementing  the  Gold  codes  as  phase  dif¬ 
fusers  for  multiplexing  a  number  of  transfer  functions  in 
a  single  hologram  are  given.  The  computations  are  carried 
out  for  Gold  codes  of  lengths  127  and  511  bits  under  the 
following  different  conditions:  The  Gold  codes  used  as  (a) 
an  ideal  phase  diffuser  with  180°  phase  difference  between 
the  elements,  (b)  an  amplitude  diffuser  with  transmittance 
values  of  0's  and  l's,  (c)  a  non-perfect  phase  diffuser  with 
phase  difference  between  the  elements  not  equal  to  180° 
and  (d)  an  ideal  phase  diffuser  illuminated  by  a  spherical 
wave  front  instead  of  a  plane  wave  front. 

The  auto  correlation  of  a  function  r^  may  be  calcu¬ 


lated  using  the  relation 
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Table  2-2.  Set  of  nine  127  bit  Gold  codes. 
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R11  "  rl*  ri  *  FCHi-Rj)  ,  (2-6) 

where  ★  represents  correlation  and  *  represents  the  conju¬ 
gate  operator  and  r ^ are  the  Fourier  transform  pairs 
defined  as 

R1  =  F{rl}  *  (2-7) 

Similarly  the  cross- correlation  between  the  functions 
r^  and  r2  may  be  obtained  from  the  relation 

R12  "  rl  *  r2  "  F{Rl‘R2}  '  (2_8) 

where  R2  =  F{r2>  (2-9) 

If  the  functions  r1  and  r2  have  spatial  widths  of  Wj^  and 
w2  respectively  then  the  functions  R^  and  R12  have  spatial 
widths  of  w^  and  W]_2  given  by 

W1X  ^  width  of  Rn  =  2WX  and 
A 

W12  =  width  of  R12  =  Wx  +  W2 

Thus  in  computer  simulations  sufficient  allowance  must  be 
made  to  accommodate  the  larger  size  of  the  output. 

Similarly  the  output  of  a  space  variant  processor 
implementing  the  phase  coded  reference  beam  technique  for 
multiplexing  may  be  simulated  by  computing  the  terms  in 
the  Eqn.  (1-14) .  Again  if  the  impulse  response  h^  has  a 


(2-10) 
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spatial  width  of  W^,  the  width  of  the  term  h^  *  r2  is 
given  by 

WQ  A  width  of  the  term  h^  *  r^  r2  *  wh  +  W1  +  W2 

(2-11) 

Thus  in  the  simulation  of  multiplex  holography  using 
Eqn.  (1-14)  the  size  of  the  output  array  must  be  made  suf¬ 
ficiently  large  as  given  by  the  above  equation. 

2.2.1  Results  of  Evaluation  of  Gold  Codes  as  Perfect  Phase 

Diffusers:  A  program  for  the  calculation  of  the  auto¬ 
correlation  and  the  cross -correlations  of  a  set  of  nine  127 
bit  Gold  codes  is  given  in  Appendix  B.  Another  program  to 
simulate  the  output  of  a  processor  using  these  codes  as 
phase  diffusers  for  multiplexing  is  given  in  Appendix  C. 

In  these  programs  the  size  of  the  output  array  is  taken  to 
be  128  elements.  Thus  in  order  to  satisfy  the  Eqn.  (2-11) 
the  widths  of  the  codes  r^  and  r2  and  the  width  of  the  im¬ 
pulse  responses  are  all  made  equal  to  42  bits.  Thus  in  these 
programs  only  the  central  42  bits  out  of  the  127  bit  codes 
are  used  in  the  computations.  (Although  it  was  possible  to 
use  64  bits  in  the  program  SPACEVAR  of  Appendix  B,  only  42 
central  bits  are  used  so  that  the  results  of  the  two  pro¬ 
grams  may  be  compared.)  The  program  SPACEVAR  computes  the 
autocorrelation  of  the  central  42  bits  of  a  127  bit  code 
and  the  cross-correlations  of  these  42  bits  with  the  central 
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42  bits  of  the  remaining  codes  in  a  set  of  9  codes,  using 
the  Eqns.  (2-6)  and  (2-8).  The  magnitudes  of  the  outputs 
are  normalized  with  reference  to  the  peak  of  the  auto-cor¬ 
relation.  The  outputs  are  plotted  to  a  width  of  2.56  inches 
and  a  height  of  2.5  inches.  The  resulting  plot  is  shown  in 
Fig.  (2-1). 

Figure  (2-1) a  is  the  plot  of  auto  correlation  of  mask 
1  and  the  Figs.  (2-1) b  through  i  are  the  cross- correlations 
of  mask  1  with  the  masks  2  through  9.  Note  that  the  cross 
correlations  have  comparable  large  magnitudes  and  hence  we 
may  expect  poor  reconstruction  of  the  impulse  responses  when 
these  codes  are  used  for  multiplexing  the  system  transfer 
functions.  The  program  MPXHOLO  of  Appendix  C  simulates  the 
output  of  a  system  when  two  transfer  functions  are  multi¬ 
plexed  using  the  Gold  codes  as  phase  diffusers.  The  pro¬ 
gram  reads  two  impulse  responses  representing  a  space-variant 
system  and  multiplexes  their  transfer  functions  on  a  single 
composite  array  using  a  different  Gold  code  as  phase  encoder 
for  each  of  the  responses.  Thus  a  composite  transfer  func¬ 
tion  hologram  is  generated.  This  hologram  includes  only 
the  terms  that  result  in  an  output  in  the  output  plane  as 
explained  in  Section  1.2.  The  output  of  the  processor  when 
the  composite  hologram  is  accessed  by  a  reference  beam 
encoded  by  a  duplicate  of  the  code  used  for  recording  is 
simulated  by  this  program.  The  program  also  simulates  the 
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output  when  only  one  of  the  transfer  functions  is  recorded 
and  played  back  using  a  Gold  code.  This  simulation  is  done 
to  assess  the  distortion  in  the  impulse  responses  due  to 
non  ideal  auto  correlation  of  the  codes.  The  program  plots 
the  outputs  of  these  simulations  as  well  as  the  impulse 
responses  used.  The  results  of  these  simulations  for 
three  sets  of  impulse  responses  are  shown  in  Figs.  (2-2) , 

(2-3)  and  (2-4) .  The  impulse  responses  used  in  the  first 
simulation  consists  of  two  disjoint  inputs  as  shown  in  Figs. 
(2-2) a  and  b.  The  outputs  when  the  transfer  functions  of 
these  impulse  responses  are  recorded  and  played  back  using 
a  Gold  code  one  at  a  time  are  shown  in  Figs.  (2-2) c  and 
d.  Note  that  due  to  non- ideal  auto-correlation  of  the  Gold 
codes,  the  output  impulse  responses  are  considerably  dis¬ 
torted.  This  distortion  is  due  solely  to  the  auto-corre¬ 
lation  and  the  effects  of  cross  correlation  are  not  included. 
Figures  (2-2)  e  and  (2-2) f  show  the  outputs  when  both  the 
transfer  functions  are  multiplexed  in  a  single  array  and 
then  an  attempt  to  retrieve  the  impulse  responses  individually 
are  made.  These  outputs  are  distorted  much  more  than  the 
outputs  of  Figs.  (2-2) c  and  d  because  of  the  cross  talk  be¬ 
tween  the  holograms  due  to  non  zero  cross- correlations  in 
addition  to  the  non  ideal  auto-correlations  of  the  Gold 
codes.  Finally  the  result  when  both  the  transfer  functions 
are  accessed  simultaneously  by  two  phase  coded  reference 
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42  central  bits  of  127  bit  Gold  e 
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beams  is  shown  in  Fig.  (2-2)g.  The  corresponding  outputs 
when  two  overlapping  impulse  responses,  one  with  a  positive 
amplitude  and  the  other  with  a  negative  amplitude,  as  shown 
in  Figs.  (2-3) a  and  b,  are  used  as  inputs  to  the  simulator 
are  shown  in  the  Figs.  (2-3) c  through  g.  In  the  next  simu¬ 
lation  two  impulse  responses  with  values  similar  to  delta 
functions  are  used  as  inputs  and  the  corresponding  outputs 
are  shown  in  Fig.  (2-4)  .  Note  that  in  this  case  the  out¬ 
puts,  although  containing  a  substantial  number  of  noise 
terms,  have  a  term  that  may  be  attributed  to  the  desired 
output.  These  simulations  show  that  the  undesired  terms 
in  the  correlations  of  the  42  central  bits  of  the  127  bit 
codes  are  too  large  and  hence  result  in  a  poor  playback. 

This  is  because  an  arbitrary  segment  of  the  Gold  code  does 
not  exhibit  the  same  degree  of  randomness  as  the  full  code. 
Thus  further  simulations  to  study  the  correlation  properties 
of  the  127  bit  codes  when  the  entire  length  of  the  codes 
are  used  for  multiplexing  were  carried  out.  This  was  done 
by  altering  the  size  of  the  arrays  in  the  computer  program 
to  accommodate  the  larger  size  of  the  outputs  as  determined 
by  the  Eqns.  (2-10)  and  (2-11).  The  results  of  these  simu¬ 
lations  are  shown  in  Figs.  (2-5)  through  (2-8).  Figure  (2-5) 
shows  the  auto  and  cross-correlations  of  the  127  bit  Gold 
code.  A  comparison  of  this  output  with  that  of  Fig.  (2-1) 
indicates  a  substantial  reduction  in  the  magnitudes  of  the 
undesired  components.  The  impulse  responses  used  in  the 
simulation  of  multiplex  holography  using  all  the  127  bits 


Figure  2-5.  Autocorrelation  and  crosscorrelations  of  127  bit  Gold 
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of  the  Gold  codes  are  similar  to  the  ones  used  in  the  simu¬ 


lation  when  only  a  part  of  the  codes  was  used  in  the  compu¬ 
tation.  Note  that  the  outputs  when  the  impulse  responses 
are  delta  like  functions  have  noise  terms  at  much  lower 
magnitudes  compared  to  the  outputs  in  Fig.  (2-4)  .  Finally 
the  results  of  evaluation  of  correlations  and  simulation  of 
multiplex  holography  when  Gold  codes  of  length  511  bits 
were  used  are  shown  in  Figs.  (2-9)  through  (2-12).  The 
Gold  codes  used  in  these  computations  are  the  outputs  of  the 
program  CODE  described  in  Section  2.1.  Here  again  note  that 
the  magnitudes  of  the  undesired  terms  in  the  correlation 
outputs  are  much  smaller  than  the  outputs  for  codes  of  smaller 
lengths.  There  is  also  a  substantial  improvement  in  the  out¬ 
put  when  delta  function  like  impulse  responses  are  used  in 
the  simulation  of  multiplex  holography  compared  to  the  out¬ 
puts  using  smaller  length  codes.  However  there  is  no  im¬ 
provement  in  the  outputs  when  the  impulse  responses  are 
broader.  This  may  be  attributed  to  the  fact  that  although 
the  magnitudes  of  the  individual  noise  elements  in  the 
correlation  outputs  are  small,  the  number  of  such  terms  are 
large  with  larger  length  of  codes  and  their  collective 
contributions  when  convolved  with  the  impulse  responses  may 
be  quite  high.  A  major  problem  may  be  the  fact  that  the  Gold 
codes,  unlike  the  maximal  length  cyclic  codes  from  which 
they  are  derived,  are  not  balanced  to  have  the  same  number 
of  +1  * s  and  -l's.  Thus  the  expected  value  of  a  bit  is  not 
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zero,  and  there  is  a  built-in  bias.  The  outputs  of  these 
simulations  suggest  that  the  method  of  space-variant 
system  representation  using  phase  coded  reference  beams  is 
more  suitable  for  space-variant  systems  having  narrow  im¬ 
pulse  responses.  An  example  of  such  a  system  is  a  magnifier 
which  transforms  points  in  the  input  plane  to  points  in  the 
output  plane. 

All  the  simulations  described  so  far  have  been  done 
under  the  assumption  that  the  Gold  code  masks  used  in  the 
system  for  recording  and  playback  are  pure  phase  masks  with 
exactly  180°  phase  difference  between  the  elements  in  the 
code.  However  it  is  generally  difficult  to  fabricate  such 
a  perfect  phase  mask.  For  this  reason  an  evaluation  of  the 
performance  of  a  system  using  non  perfect  phase  masks  and 
binary  amplitude  masks  was  carried  out  and  the  results  are 
presented  in  the  next  subsections . 

2.2.2  Results  of  Evaluation  of  Gold  Codes  as  Amplitude 
Masks  and  Non  Perfect  Phase  Masks:  An  amplitude 
mask  of  a  binary  Gold  code  may  be  easily  fabricated  using 
any  high  contrast  copy  film.  The  transmittance  of  these 
masks  has  values  of  0  and  1  instead  of  +1  and  -1  for  a  per¬ 
fect  phase  mask.  The  auto-correlation  and  the  cross-corre¬ 
lations  of  such  amplitude  masks  with  code  lengths  of  127 
bits  are  shown  in  Fig.  (2-13).  Note  that  the  magnitudes 
of  the  correlation  terms  are  very  large  toward  the  center 
of  the  output  array  and  then  tail  off  rather  slowly.  The 
outputs  of  simulation  of  multiplex  holography  for  the  same 
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inputs  used  in  the  case  of  pure  phase  masks  are  shown  in 
Figs.  (2-14)  through  (2-16) .  Comparison  of  these  outputs 
with  those  of  Figs.  (2-5)  through  (2-8)  reveals  that  the 
phase  masks  are  much  superior  in  performance  to  amplitude 
masks.  However  as  mentioned  earlier  it  is  difficult  to  fab¬ 
ricate  phase  masks  with  phase  difference  of  exactly  180 
degrees  between  the  elements.  Thus  a  study  was  made  to 
determine  an  acceptable  level  of  tolerance  in  the  value 
of  the  phase  difference.  The  program  SPACEVAR  of  Appendix 
B  was  modified  to  account  for  non  perfect  phase  masks  with 
phase  differences  of  172°,  162°,  150°,  and  120°.  The 
output  of  auto- correlation  of  a  127  bit  mask  with  itself  and 
the  cross-correlation  with  two  other  127  bit  masks  were 
computed  and  plotted.  The  plots  for  various  degrees  of  non 
perfectness  is  shown  in  Fig.  (2-17) .  From  these  outputs 
it  may  be  concluded  that  it  is  desirable  that  the  phase 
differences  between  the  elements  of  the  code  should  be 
within  10%  of  180°.  In  the  next  subsection  the  results 
of  evaluation  of  Gold  codes  when  a  spherical  wavefront  is 
used  during  recording  and  playback  are  presented. 

2.2.3  Results  of  Evaluation  of  Gold  Codes  Illuminated  by 

Spherical  Wavefront;  It  has  been  reported  [16]  that 
the  use  of  spherical  wave  illumination  instead  of  plane 
wave  illumination  in  an  optical  processing  system  imple¬ 
menting  the  phase  coded  reference  beam  multiplexing  tech- 
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nique  results  in  a  reduction  in  the  magnitudes  of  the  cross¬ 
talk  terms  during  playback.  A  computer  program  SPHWAVE 
given  in  the  Appendix  D  was  used  to  compute  the  auto  and 
cross-correlations  of  two  127  bit  masks  when  illuminated 
by  spherical  waves  of  different  values  of  chirp,  i.e., 
different  radii  of  curvature  and  the  width  of  mask. 

The  method  for  calculating  the  chirp  at  each  element 
of  the  mask  is  shown  in  the  Fig.  (2-18) .  Let  R  be  the  radius 
of  curvature  of  spherical  wave  in  millimeters,  and  let  w 
be  the  width  of  the  mask  in  millimeters;  the  path  difference 
M  between  the  wave  front  at  a  point  n  elements  from  the 
center  of  the  array  is  then  given  by 

M  a  /r2  +(!gW)2  -  R  ,  (2-12) 

where  N  is  the  total  number  of  elements  in  the  entire  array. 

The  phase  difference  0  in  radians  at  the  center  of 
the  element  with  reference  to  the  center  of  the  array  when 
using  an  illumination  of  wavelength  equal  to  X  millimeters 
is  given  by 

0  =  (i*  -  n)  2tt  ,  (2-13) 

where  n  is  an  integer  chosen  such  that  0  _<  0  <  2ir. 

At  optical  wavelengths  the  phase  angle  8  changes  very 
rapidly  along  the  width  of  the  phase  mask.  There  will  be 
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many  complete  cycles  of  phase  change  within  each  element 
of  a  127  bit  mask  even  for  such  small  mask  sizes  as  3  mms. 
For  this  reason  it  is  necessary  to  breakup  each  element 
of  the  code  into  several  subcells  as  is  done  in  the  pro¬ 
gram.  Then  the  value  of  each  subcell  is  determined  by  the 
value  of  the  originating  element,  corrected  for  the  phase 
change  due  to  the  shperical  wave  front  at  the  center  of  the 
subcell.  The  auto-correlation  output  of  a  127  bit  code  and 
its  cross-correlations  with  another  127  bit  code  for  differ¬ 
ent  values  of  chirp  as  determined  by  the  radius  of  wavefront 
and  the  width  of  mask  are  shown  in  Figs.  (2-19) a  and  (2-19) b 
respectively.  These  plots  are  scaled  in  width  to  account 
for  the  variation  in  the  size  of  the  masks.  The  heights 
are  scaled  so  that  the  areas  under  each  of  the  auto  corre¬ 
lation  peaks  are  equal  in  order  to  establish  a  criterion 
for  comparison.  Note  that  the  illumination  by  a  spherical 
wavefront  has  a  tendency  to  reduce  the  amount  of  undesired 
terms  in  the  correlation  that  are  located  away  from  the  cen¬ 
ter  of  the  peak  of  the  auto- correlation .  The  terms  near  the 
center  are  not  changed  appreciably.  Also  note  that  as  the 
radius  of  wavefront  gets  very  large  the  correlation  outputs 
approach  that  of  a  mask  illuminated  by  a  plane  wave  as  ob¬ 
tained  in  Fig.  (2-1) .  Further  simulations  are  necessary 
to  quantitatively  establish  the  exact  amount  of  improvement 
in  the  output  that  may  be  obtained  by  using  the  spherical 


Figure  2-19.  Autocorrelation  and  crosscorrelations  of  42 

central  bits  of  127  bit  Gold  codes  illuminated 
by  a  spherical  wave  of  radius  R  and  width  of 
mask  W. 
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wave  front  as  well  as  for  determining  the  optimum  relations 
between  the  width  of  the  mask  and  the  radius  of  curvature 
of  the  spherical  wave  front. 

The  results  of  the  computer  simulations  presented 
in  this  chapter  are  useful  in  choosing  a  space-variant 
system  for  representation  using  the  phase  coded  reference 
beam  approach.  The  need  for  a  good  phase  mask  with  phase 
differences  between  the  elements  close  to  180°  was  also 
established.  A  technique  for  fabricating  a  two  dimensional 
phase  mask  using  Dichromated  gelatin  is  described  in 
Appendix  G. 

In  the  following  two  chapters  an  alternate  multiplex¬ 
ing  technique  for  generating  a  composite  transfer  function 
hologram  is  presented 


CHAPTER  3 


SAMPLED  INPUT/ SAMPLED  TRANSFER  FUNCTION  APPROACH 

3.1  Space  Division  Multiplexing  of  Transfer  Functions 

Consider  a  system  sampled  at  N  points  in  the  input 
plane,  as  determined  by  the  sampling  theorem  of  Eqn.  (1-10) 
As  a  result  we  have  N  line  spread  functions  represent¬ 
ing  the  system.  Thus  after  Fourier  transformation  there 
are  N  transfer  functions  in  the  holographic  plane  to  be 
multiplexed  in  a  single  recording  medium.  When  the  system 
line  spread  functions  are  space  limited,  it  is  possible  to 
sample  their  transfer  functions  at  a  rate  determined  by 
the  modified  version  of  the  Whittaker-Shannon  sampling 
theorem  (11]  and  generate  a  composite  hologram  containing 
the  samples  of  all  the  transfer  functions.  A  typical 
space  limited  line  spread  function  might  be  as  shown  in 
Fig.  (3-1)  .  The  maximum  spatial  width  of  this  function 
is  2x*  where  x*  is  the  larger  of  the  values  on  either  side 

mm 

of  the  axis  of  the  optical  system.  Let  XM  be  the  maximum 
of  {x^}  for  all  i  *  1,2,...N.  Then  the  maximum  sampling 
interval  required  in  the  transfer  function  plane  is  given 
by  the  Whittaker-Shannon  sampling  theorem: 


where  Afx  has  the  dimensions  of  spatial  frequency;  or  in 
terms  of  linear  dimensions, 
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AU  *  AfAf 


(3-2) 


Where  A  is  the  wave  length  of  coherent  light,  f  is  the  focal 
length  of  the  Fourier  transforming  lens  and  AU  is  the  samp¬ 
ling  interval  in  the  holographic  plane  in  units  of  length. 
Since  we  have  N  transfer  functions  to  be  multiplexed,  the 
width  of  each  individual  sample  is  given  by 


Av  ■  M.  .  (3-3) 

N 

Again  if  the  magnitudes  of  the  transfer  functions  are 
essentially  zero  beyond  the  width  |u|  _>  b/2,  the  transfer 
functions  may  be  approximately  represented  by  limiting  the 
number  of  samples  to 


m  ,  b,  s 

AU  Af 


(3-4) 


instead  of  the  infinite  number  of  samples  required  by  the 
sampling  theorem. 

An  example  of  sampling  a  transfer  function  and  the 
spatial  distribution  of  samples  in  the  holographic  plane, 
for  N  =  3  and  M  =  4,  is  shown  in  Fig.  (3-2)  .  Each  sample 
is  marked  as  Hj,  when  i  represents  the  ith  transfer  func¬ 
tion  being  multiplexed  and  j  represents  the  jth  sample  of 
the  itJl  transfer  function. 

3.2  An  Optical  Recording  And  Playback  Scheme 

A  scheme  for  implementing  the  multiplexing  technique 
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just  described  is  shown  in  Fig.  (3-3)  .  During  recording, 
a  binary  mask  with  the  width  of  each  of  the  N  transparent 
areas  equal  to  Av  spaced  at  intervals  of  AU  is  placed 
immediately  in  front  of  the  recording  medium.  This  mask 
samples  the  transfer  function  hologram  at  intervals  of 
AU.  The  mask  is  moved  by  a  distance  of  Av  after  recording 
each  hologram.  Thus  at  the  end  of  recording  and  process¬ 
ing,  assuming  that  the  resultant  transmittance  after  pro¬ 
cessing  is  proportional  to  intensity,  the  transmittance 
of  the  hologram  is  given  by 


N 

t  (u)  =  l  [|H1(u)  +  R(u)  |  2]  [Rect  (“_)  *  Comb  ( -^A-)  1 

(3-5) 


where  *  represents  convolution, 


Comb(x)  =  l  6 (x-n)  (3-6) 

n=-<» 

and  Rect  (x)  is  as  defined  in  equation  (1-8)  . 

For  playing  back  this  multiplexed  hologram  the  scheme 
is  as  shown  in  Fig.  (3-4)  .  The  binary  mask  used  in  the 
recording  step  is  now  replaced  by  a  multiple  input  function 
transparency  with  transmittance  equal  to  f(C^)  at  all  M 
points  that  are  directly  in  front  of  the  samples  of  the 
ith  transfer  function.  The  transmittance  of  this  trans¬ 
parency  may  be  represented  as 

N 

s(u)  =  l  f  (Ci)Rect(H_)  *  Comb  (U.~AA.V.)  .  (3-7) 

1  Av  AU 


When  this  transparency  is  illuminated  by  the  reference 


beam  R,  the  reconstructed  wavefront  to  the  right  of  the 
hologram  is  given  by 


G  (u)  =  R  (u)  S  (u)  t  (u) 


N 

«  £^R(  |  H1  +  R|  2)  f  (Ci)Rect(~)  *  Comb  . 

(3-8) 

Now  out  of  the  four  terms  in  the  expansion  of  r'h^  +  r|2 

* 

(Eqn.  (1-13))  only  the  term  RR  is  diffracted  by  the  holo¬ 
gram  in  the  direction  of  the  output  plane  as  a  result  of 
the  offset  in  the  reference  beam.  Thus  the  output  after 
Fourier  transformation  by  the  lens  L2  is  given  by 


g’  (x) 


N 


F"1[.?iHiRR*f  (5i)Rect(^?)  *  Comb  (Hl^)  ] 


N 


l  K  h.  *  (r  *  r)  *  [f  (5.)  Sine  (|_)  (Comb  (^f-)  ] 
i=l  1  1  ^nxM  ^XM 


(3-9) 


where  ★  represents  correlation. 

Here  K  is  a  scaling  factor  due  to  Fourier  transforma¬ 
tion. 

In  this  equation  the  term  r  ★  r  approaches  a  delta  func¬ 
tion  if  the  reference  source  r  approaches  a  delta  function. 
The  term  Sine  ( 2~*-x-)  ^ue  to  the  ^inite  size  of  the  sample 
width  in  the  frequency  plane,  and  approaches  a  constant  in 
the  limit  as  Av->-0.  Finally  the  term  Comb  )  is  present 
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because  of  the  sampling  carried  out  for  multiplexing  the 
holograms.  Thus  under  the  assumption  that  the  width  of  each 
sample  Av  is  small,  and  that  r  is  a  delta  function,  the 
equation  for  g'  (x)  is  given  by 

N 

g'(x)  =  l  f(5<)h.  *  Comb (•=£—)  (3-10) 

i=l  1  1  ^XM 

The  required  output  g(x)  is 

N 

g(x)  =  I  f(5,)h.  (3-11) 

i=l  1  1 

Hence  to  recover  the  output  g(x)  from  g'  (x)  we  need  a 
mask  in  the  output  plane  with  a  slit  which  passes  one  of 
the  multiple  images. 

All  the  mathematical  derivations  carried  out  so  far 
has  been  in  1-D  for  clarity  of  presentation.  Extensions 
to  two  dimensions  are  straightforward. 

This  method  of  multiplexing  does  not  require  the 
multiple  reference  beams  that  were  necessary  in  the  encoded 
reference  beam  approach  described  in  Section  1.2.  However 
the  need  exists  for  preparing  a  multiple  input  function  mask 
for  the  spacelimited  input  f(£).  This  mask  is  used  during 
the  playback  as  explained  in  the  previous  paragraphs.  This 
mask  generates  coherent  replications  of  the  input  function 
to  illuminate  the  hologram.  Some  of  the  schemes  for 
achieving  this  objective  are  described  in  Section  3.4.  In 
the  next  section  the  results  of  1-D  computer  simulations 


carried  out  to  verify  the  sampled  transfer  function  multiplex¬ 
ing  technique  are  described. 


3.3  One-Dimensional  Computer  Simulations 

The  computer  program  used  to  simulate  the  sampled 
transfer  function  multiplexing  technique  is  given  in 
Appendix  E.  In  this  simulation  the  number  of  samples  in 
the  input  plane  is  taken  as  N  =  4.  Thus  there  are  four 
transfer  functions  to  be  multiplexed  in  a  single  hologram. 

A  one  dimensional  array  of  128  elements  is  used  to  repre¬ 
sent  each  impulse  response.  This  array  is  used  as  the  in¬ 
put  to  the  discrete  Fast  Fourier  Transform  routine  (FFT) 
to  generate  a  128  element  array  of  Fourier  components.  This 
transfer  function  array  is  sampled  at  an  interval  of  four 
elements  and  the  samples  are  stored  at  their  corresponding 
positions  in  another  array  representing  the  composite  holo¬ 
gram.  The  above  steps  are  repeated  for  all  the  four  impulse 
responses,  resulting  in  a  final  composite  array  of  128  ele¬ 
ments.  This  composite  array  is  played  back  by  using  it  as 
the  input  to  a  second  Fourier  transform  routine.  As  the 
composite  array  is  not  multiplied  by  any  term  representing 
the  input  function,  the  result  after  Fourier  transformation 
should  be  the  sum  of  the  individual  impulse  responses.  In 
fact  this  simulation  is  equivalent  to  a  situation  when  the 
input  is  a  constant  for  all  the  impulse  responses.  The 
output  of  the  system  when  the  transfer  functions  are  accessed 
one  at  a  time  is  also  simulated. 


59 


The  four  impulse  responses  used  in  this  experiment  are 
shown  in  Fig.  (3-5a) .  Note  that  since  we  have  four  trans¬ 
fer  functions  to  be  multiplexed  in  a  composite  array  of  128 
elements,  the  maximum  extent  xM  of  any  of  these  impulse  re¬ 
sponses  should  be  less  than  16  elements  on  either  side  of 
the  center  of  the  array,  i.e.,  2xM  <  128/4  *  32  elements; 
otherwise  aliasing  problems  will  result  in  the  output  plane 
when  the  multiplexed  transfer  function  array  is  played  back. 

The  magnitudes  of  the  transfer  functions  of  each  of  these 
impulse  responses  are  shown  in  Fig.  ( 3— 5b) .  These  transfer 
functions  are  sampled  by  selecting  the  1st,  5^,  9^  .... 
elements  of  the  first  transfer  function,  2nd,  6th  ....  ele¬ 
ments  of  the  second  transfer  function,  3rd,  7th  ....  elements 
of  the  third  transfer  function  and  4th,  8th  ....elements  of 
the  last  transfer  function.  These  samples  are  placed  in 
their  respective  positions  in  another  composite  array.  The 
magnitude  of  the  elements  in  this  array  is  shown  in  Fig. 

(3-6a) .  This  array  is  then  Fourier  transformed  and  the  out¬ 
put  is  shown  in  Fig.  (3-6b) .  This  output  represents  the 
system  output  when  illuminated  by  an  input  function  f(£) 
which  is  a  constant  over  all  the  sample  points  in  the  in¬ 
put  plane.  Note  that  the  output,  which  is  a  sum  of  all  the 
four  impulses  responses,  is  replicated  four  times  in  the 
output  plane.  This  is  because  of  the  Comb  (x/2xM)  term  in 
Eqn.  (3-10).  The  coordinate  reversal  is  due  to  the  consecu¬ 
tive  application  of  two  Fourier  transformation  operations.  Next, 
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to  simulate  individual  access  of  each  impulse  response,  the 
samples  belonging  to  each  transfer  function  are  recovered 
from  the  composite  array  and  are  placed  in  another  array. 
The  plot  of  the  magnitudes  of  the  elements  in  this  array 
is  shown  in  Fig.  (3-7a) .  The  results  of  Fourier  trans¬ 
forming  these  arrays  are  shown  in  Fig.  (3-7b) .  This  simu¬ 
lation  is  equivalent  to  the  situation  when  the  input  func¬ 
tion  f (C)  has  a  magnitude  of  1  at  the  input  sample  point 
corresponding  to  the  impulse  response  being  accessed  and 
zero  at  all  the  other  sample  points.  Again  note  that  the 
output  is  replicated  four  times  in  the  output  plane.  This 
program  was  written  in  the  Fortran  63  language  for  use  with 
a  CDC  1604  Computer  and  Cal-Comp  drum  type  plotter. 

3.4  Schemes  for  the  Generation  of  Multiples  of  the  Input 
Function 

In  this  section  some  of  the  schemes  for  generating 
multiple  images  of  the  input  function  are  briefly  described 

3.4.1  Multiple  Image  Transparencies;  Multiple  copies  of 
the  sampled  values  of  the  input  function  f (£)  are  prepared 
on  a  film  using  a  step-and-repeat  process  or  other  methods. 
A  typical  transparency  for  a  given  function  f(C)  is  shown 
in  Fig.  (3-8).  The  disadvantages  of  this  method  are  (a) 
Separate  masks  are  required  for  each  of  the  input  functions 
(b)  the  process  is  slow. 
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3.4.2  Multiple  Imaging  Using  Beam  Splitters:  An  array  of 
beam  splitters  may  be  set  up  to  generate  multiple  images 
of  the  input  function.  The  requirement,  however  that  all 
the  images  should  have  the  same  amplitude  and  be  coherent 
with  one  another  requires  high  precision  in  the  values 

of  transmittance  and  reflectance  of  each  of  the  elements 
as  well  as  in  the  optical  path  lengths. 

3.4.3  Multiple  Imaging  Using  Phase  Holograms:  The  use  of 
phase  holograms  to  produce  equally  bright  multiple  images 
in  the  fabrication  of  integrated  circuits  has  been  reported 
[17] .  Similar  techniques  to  produce  multiple  coherent  images 
may  be  possible. 

3.4.4  Use  of  Fiber  Optic  Elements:  Bundles  of  equal  lengths 
of  fiber  optic  elements  may  be  arranged  as  shown  in  Fig. 

(3-9)  to  generate  multiple  images  of  the  sampled  input  func¬ 
tion  when  illuminated  by  a  plane  wave. 

3.4.5  Use  of  Liquid  Crystal  Devices;  The  property  of  liquid 
crystal  devices  by  which  local  changes  in  the  periodicity 

of  a  phase  grating  becomes  proportional  to  the  light 
variation  incident  on  the  device  has  been  used  in  optical 
computing  [18].  This  property  may  be  used  in  conjunction 
with  other  techniques  described  earlier  to  generate  coherent 
multiple  images. 
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In  the  next  chapter  the  experimental  results  obtained 
using  2-D  computer  multiplexed  holograms  are  presented. 


CHAPTER  4 


EXPERIMENTAL  RESULTS  USING 
COMPUTER  MULTIPLEXED  HOLOGRAMS 

The  technique  of  representing  a  space-variant  system 
using  a  thin  recording  medium  has  the  advantage  of  allow¬ 
ing  computer  generation  of  the  system  transfer  function 
hologram  instead  of  using  the  optical  recording  schemes 
described  in  section  (3-2) .  In  this  chapter  the  results  of 
2  experiments  using  computer  multiplexed  holograms  are  pre¬ 
sented.  The  experiments  were  conducted  to  verify  the  tech¬ 
nique  of  multiplexing  the  transfer  functions  using  the  samp¬ 
ling  method  described  in  Chapter  3.  In  these  experiments 
the  number  of  samples  N  in  the  input  plane  was  taken  to  be 
2  in  each  dimension,  thus  requiring  the  multiplexing  of  N  x  N 
4  transfer  functions  on  a  single  hologram. 

4.1  Computer  Generation  and  Playback  of  the  Multiplexed 
Hologram 

The  computer  program  used  to  generate  the  transfer 
function  hologram  and  to  simulate  the  playback  of  the  system 
is  given  in  Appendix  F.  In  this  program  a  2-D  array  of 
64  x  64  elements  is  used  to  represent  an  impulse  response. 

The  four  impulse  responses  are  Fourier  transformed  using 
the  discrete  Fast  Fourier  Transform  (FFT)  subroutine  to 
result  in  four  transfer  function  arrays  of  64  x  64  elements 
each.  These  four  transfer  function  arrays  are  multiplexed 
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into  a  single  composite  array  of  64  x  64  elements  by  select¬ 
ing  every  alternate  element  from  each  of  the  arrays  in  both 
direction  as  shown  in  Fig.  (4-1)  i.e.,  the  1st,  3rd,  5th, 
etc.  element  from  the  1st,  3rd,  5&1,  etc.,  rows  of  the  trans¬ 
fer  function  are  selected  and  placed  in  their  correspond¬ 
ing  positions  in  the  composite  array.  Similarly  the  2nd, 

4fch,  6th,  etc.  element  from  the  1st,  3rd,  5th,  etc.  rows 
of  the  transfer  function  H 2,  1st,  3rd,  5th,  etc.  elements 
from  the  2nd,  4^,  6^,  etc.  rows  of  the  transfer  function 
H3  and  2nd,  4^,  6t*l#  etc.  elements  from  the  2nd,  4^,  6^h/ 
etc.  rows  of  the  transfer  function  H4  are  selected  and  are 
placed  in  their  corresponding  positions  in  the  composite 
array.  Thus  in  this  scheme  the  sampling  interval  in  the 
transfer  function  plane  is  2  elements.  This  implies  that 
the  maximum  size  of  the  impulse  response  in  either  direc¬ 
tion  from  the  center  of  the  array  must  be  less  than  64/ (2x2) 

=  16  elements  in  order  to  satisfy  the  sampling  theorem. 

In  general  if  N  transfer  functions  are  to  be  multiplexed 
in  each  dimension  the  maximum  extent  of  the  impulse  re¬ 
sponse  from  the  center  of  the  array  is  given  by 

y. 

XM  =  2N  elements,  (4-1) 

where  b  is  the  number  of  elements  in  the  array  in  each 
dimension.  If  this  condition  is  not  satisfied  aliasing 
errors  will  result  during  the  playback  step. 


(a)  Selection  of  samples  from  the  transfer  functions. 


(b) 


Figure 


Position  of  components  from  the  fc  r  transfer 
functions  in  the  composite  array. 


4-1.  Scheme  for  computer,  multiplexing  the  transfer 
functions . 
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At  the  end  of  the  multiplexing  step  we  have  one  com¬ 
posite  array  of  64  x  64  elements  representing  the  multi¬ 
plexed  hologram.  The  magnitude  and  phase  of  the  elements 
in  this  array  are  plotted  using  Burckhardt's  3  vector  method 
[19]  in  cells  of  size  0.15"  x  0.15",  resulting  in  a  plot  of 
size  9.6"  x  9.6".  Since  the  resolution  of  the  plotter  was 
limited  to  0.01"  the  magnitudes  of  the  elements  in  the  array 
are  quantized  to  a  total  of  15  levels.  Thus  all  elements 
whose  magnitudes  are  less  than  l/15th  of  the  magnitude  of 
the  largest  element  in  the  array  are  set  to  zero.  These 
quantized  vectors  are  resolved  into  components  along  3 
vectors  120°  apart  as  shown  in  Fig.  (4-2) .  Each  of  these 
components  are  represented  on  the  plot  by  three  subcells 
of  0.05"  width.  The  height  of  these  subcells  is  propor¬ 
tional  to  the  magnitude  of  the  component.  The  plot  of  one 
such  cell  corresponding  to  the  element  having  magnitude 
and  phase  as  shown  in  Fig.  (4-2)  is  given  in  Fig.  (4-3) . 

The  plot  of  the  encoded  sampled  transform  array  is 
reduced  to  a  size  of  0.4"  x  0.4"  using  high  contrast  copy 
film  and  is  used  in  the  optical  system  shown  in  Fig.  (4-4)  . 
The  optical  playback  system  consists  of  lens  placed  at 
a  distance  of  one  focal  length  from  both  the  hologram  and 
the  output  plane.  The  hologram  is  placed  in  the  plane  U 
and  the  output,  after  Fourier  transformation  by  lens  , 
appears  in  the  X  plane.  A  photograph  of  the  optical  setup 
is  given  in  Fig.  (4-5)  .  In  this  setup  an  additional  lens  is 
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Figure  4-5.  Photograph  of  the  Optical  System  used  for  the 
playback  of  computer  multiplexed  holograms. 
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used  to  project  enlarged  images  in  the  output  plane.  The 
entire  hologram  is  illuminated  by  a  plane  wave  of  constant 
amplitude.  This  is  equivalent  to  the  situation  when  the 
input  function  ftC^  in  Eqn.  (3-10)  is  constant  for  all  the 
impulse  responses  hj_,  i  *  1,2, 3, 4.  The  output  is  observed 
as  an  intensity  distribution  in  the  output  plane.  The 
computer  program  also  simulates  the  playback  of  the  optical 
system.  This  is  done  by  using  the  composite  transfer  func¬ 
tion  array  as  the  input  to  the  Fast  Fourier  Transform  routine 
and  the  magnitude  of  the  output  is  plotted  as  before. 

Next,  to  simulate  the  playback  of  only  one  impulse 
response  the  components  belonging  to  one  of  the  transfer 
function  from  the  composite  array  are  selected  and  placed 
in  their  respective  position  in  another  array  with  magni¬ 
tudes  of  all  other  elements  set  to  zero.  A  hologram  of  this 
transfer  function  array  is  prepared  as  before  using  a  high 
contrast  copy  film.  This  hologram  is  used  in  the  optical 
system  of  Fig.  (4-4)  and  the  output  is  observed  as  an  in¬ 
tensity  distribution  in  the  output  plane.  This  is  equiva¬ 
lent  to  the  situation  when  the  input  function  f  ( Ci )  in 
Eqn.  (3-10)  is  nonzero  at  only  the  point  corresponding  to 
the  impulse  response  h^  being  accessed  and  zero  at  all  other 
points.  In  the  following  sections  the  results  obtained 
using  this  computer  multiplexing  technique  are  presented. 


The  four  disjoint  functions  representing  the  impulse 
responses  used  in  this  experiment  are  shown  in  Pig.  (4-6)  . 
Note  that  the  maximum  extent  of  these  impulse  responses  in 
either  x  or  y  direction  with  respect  to  the  center  of  the 
array  is  less  than  16  elements  and  hence  satisfies  the  con¬ 
straints  imposed  by  Eqn.  (4-1)  .  These  impulse  responses  are 
used  as  inputs  to  the  computer  program  described  in  the 
previous  section.  The  composite  hologram  generated  by  the 
computer  is  shown  in  Fig.  (4-7)  .  The  result  of  playback  of 
this  hologram  in  the  optical  system  of  Fig.  (4-4)  is  shown 
in  Fig.  (4-8)  .  The  binary  mask  shown  in  Fig.  (4-4)  was  not 
used  while  recording  this  output.  The  multiple  outputs 
seen  in  this  output  are  due  to  two  reasons.  First,  as  a 
result  of  the  sampling  in  the  transfer  function  plane 
multiple  images  are  produced  in  the  output  plane  as  illus¬ 
trated  by  Eqn.  (3-10)  .  Second,  the  Fast  Fourier  Transform 
subroutine  assumes  that  the  object  at  the  input  is  one 
period  of  a  periodic  function  in  both  x  and  y  directions  so 
that  the  output  is  limited  to  the  size  of  one  period  of  the 
array.  The  result  of  optical  playback  using  the  binary 
mask  as  shown  in  Fig.  (4-4)  to  pass  only  one  of  the  multiple 
images  is  shown  in  Fig.  (4-9).  Note  that  this  output  is  a 
sum  of  all  the  four  impulse  responses.  This  simulation  is 
equivalent  to  accessing  all  the  impulse  responses  by  an 


Figure  4-7.  Computer  multiplexed  hologram  using  impulse 
responses  of  Fig.  4-6. 


Figure  4-8.  Output  of  the  Optical  System 
when  the  hologram  of  Fig.  4-7 
is  played  back. 
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Figure  4-9.  Enlarged  output  of  the  Optical  System 

when  the  hologram  of  Fig.  4-7  is  played 
using  a  binary  mask  to  pass  only  one  of 
the  multiple  images. 
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Figure  4-10.  Computer-simulated  output  when  all  the  impulse 
responses  of  Fig.  4-6  are  played  back. 


Figure  4-11.  Output  of  the  Optical  System  when 
only  one  of  the  impulse  responses 
of  Fig.  4-6  is  played  back. 


Output  of  the  Optical  System  when  one 
of  the  impulse  responses  of  Fig.  4-6 
is  played  back  using  a  binary  mask  to 
pass  only  one  of  the  multiple  images. 


Figure  4-12. 
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input  function  which  is  constant  at  all  the  sample  points 
in  the  input  plane.  The  result  of  a  computer  simulation  of 
the  playback  is  shown  in  Fig.  (4-10)  .  As  explained  in  the 
previous  section  the  components  belonging  to  one  of  the 
transfer  functions  are  selected  from  the  composite  array  and 
a  hologram  is  prepared.  The  result  of  playback  of  this  holo¬ 
gram  is  shown  in  Fig.  (4-11)  .  Again  an  enlarged  output  of 
one  of  the  multiple  images  in  the  output  plane  is  shown  in 
Fig.  (4-12) .  This  simulation  is  equivalent  to  the  situation 
when  the  input  function  accesses  only  one  of  the  impulse 
responses . 

4.3  Experimental  Results  Using  Overlapping  Impulse  Responses 
The  impulse  responses  used  in  this  experiment  are  shown 
in  Fig.  (4-13) .  The  first  and  second  functions  are  made  up 
of  two  disjoint  right  angle  triangles  with  the  values  of 
elements  within  these  triangles  equal  to  -1  and  zero  every¬ 
where  else.  Similarly  the  elements  within  the  inverted 
triangles  of  the  3rc^  and  4th  function  have  a  value  of  +1.  It 
is  clear  that  the  set  of  functions  3  and  4  partially  overlap 
the  set  of  functions  1  and  2.  Thus  when  all  the  functions 
are  added  together  a  central  hexogonal  area  of  zeros  are 
generated  surrounded  by  a  star  like  outer  pattern.  This  ex¬ 
periment  was  conducted  to  verify  whether  the  property  of 
coherent  addition  is  retained  when  the  transfer  functions 
are  multiplexed  using  the  sampled  transfer  function  approach. 
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Figure  4-14.  Computer  multiplexed  hologram  using  impulse  responses 
of  Fig.  4-13. 
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Figure  4-16.  Computer  simulated  output  when  all  the  impulse 
responses  of  Fig.  4-13  are  played  back. 
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Output  of  the  Optical  System  when 
only  one  of  the  impulse  responses 
of  Fig.  4-13  is  played  back. 
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The  plot  of  the  composite  transfer  function  is  shown  in 
Pig.  (4-14) .  The  result  of  optical  playback  of  this  com¬ 
posite  transfer  function  is  shown  in  Fig.  (4-15)  and  the 
computer  simulated  output  is  shown  in  Pig.  (4-16)  .  Note 
that  the  output  contains  a  central  hexogonal  array  of  zeros 
and  hence  verifies  the  coherent  addition  property.  The 
result  of  playback  of  one  of  the  impulse  responses  is 
shown  in  Fig.  (4-17) . 

4.4  Multiplication  of  Impulse  Responses  by  a  Phase  Function 
It  was  mentioned  in  section  4-1  that  some  of  the  terms 
in  the  composite  transfer  function  array  are  set  to  zero  if 
their  magnitude  is  less  than  1/15  times  the  magnitude  of 
the  largest  component  in  the  array.  This  results  in  the 
loss  of  many  terms  especially  when  the  transfer  function  has 
a  dominating  term  of  very  large  magnitude  (usually  the  zero 
frequency  term) .  This  results  in  a  computer  hologram  plot 
with  only  a  few  terms.  The  playback  of  such  a  hologram  re¬ 
sults  in  poor  reconstruction.  To  circumvent  this  problem 
the  impulse  responses  are  multiplied  by  a  phase  function 
having  magnitudes  of  +1  and  -1  arranged  in  a  checkerboard 
pattern  as  shown  in  Fig.  (4-18) .  In  general  this  operation 
spreads  the  components  in  the  Fourier  plane  more  evenly  and 
the  plot  of  the  hologram  after  quantization  contains  more 
components.  However  multiplication  of  the  impulse  responses 
by  this  phase  function  will  not  result  in  any  change  in  the 


Figure  4-19.  Composite  hologram  when  the  impulse  responses 
of  Fig.  4-6  are  multiplied  by  the  phase  mask 
of  Fig.  4-18. 
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Figure  4-21.  Computer  simulated  output  when  all  the  impulse 
responses  of  Fig.  4-6  premultiplied  by  the 
phase  mask  of  Fig.  4-18  are  played  back. 
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observed  distribution  in  the  output  plane  during  playback 
since  the  output  is  observed  as  an  intensity  distribution. 
Also  a  better  reconstruction  is  obtained  as  more  components 
are  present  in  the  hologram.  The  plot  of  the  composite 
transfer  function  corresponding  to  the  impulse  responses 
of  Fig.  (4-6)  using  the  phase  mask  as  premultiplier  is  shown 
in  Fig.  (4-19)  .  One  may  compare  this  plot  with  that  of  Fig. 
(4-7)  which  was  produced  without  premultiplication  by  the 
phase  mask.  The  result  of  optical  playback  of  this  holo¬ 
gram  is  shown  in  Fig.  (4-20)  .  The  result  of  computer  simu¬ 
lation  of  this  output  is  shown  in  Fig.  (4-21) .  Similarly 
the  plot  of  the  composite  transfer  function,  the  optical 
playback  of  the  composite  transfer  function  and  the  com¬ 
puter  simulated  output  when  the  impulse  responses  of  Fig. 
(4-13)  are  multiplied  by  the  phase  function  of  Fig.  (4-18) 
are  shown  in  Figures  (4-22) ,  (4-23) ,  and  (4-24)  respectively. 
The  block  like  structure  of  the  output  pattern  could  be  re¬ 
moved  by  using  a  random  phase  premultiplexing  mask  rather 
than  a  periodic  one. 

4.5  Computer  Multiplexing  Using  Low  Pass  Filtered  Transfer 
Functions 

When  the  magnitudes  of  the  high  spatial  frequency  com¬ 
ponents  of  the  transfer  functions  representing  the  space- 
variant  system  are  small,  the  quantized  composite  transfer 
function  array  contains  components  only  in  the  central  region 
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Figure  4-23.  Output  of  the  Optical  System  when 

the  hologram  of  Fig.  4-22  is  played 
back . 
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Figure  4-24.  Computer  simulated  output  when  all  the  impulse 
responses  of  Fig.  4-13  premultiplied  by  the 
phase  function  of  Fig.  4-18  are  played  back. 
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In  such  a  situation  a  different  scheme  for  selecting  the 
samples  from  the  transfer  function  arrays  for  multiplexing 
maybe  employed.  In  this  scheme  the  central  32  x  32  elements 
from  each  of  the  transfer  function  arrays  through  H4  are 
selected  and  repositioned  in  a  composite  transfer  function 
array  of  64  x  64  elements  as  shown  in  Fig.  (4-25)  .  This 
scheme  is  different  from  the  one  described  in  section  (4-1) 
in  which  alternate  elements  in  both  directions  were  selected 
as  samples  and  were  placed  in  their  corresponding  positions 
in  the  composite  array.  The  present  scheme  is  equivalent  to 
low  pass  filtering  the  transfer  functions  and  hence  infor¬ 
mation  about  high  frequency  components,  such  as  sharp  edges, 
is  lost.  Since  all  the  elements  in  the  center  of  the  trans¬ 
fer  function  array  are  used,  however,  effectively  no  sampling 
is  being  done  in  the  Fourier  plane,  and  hence  the  restriction 
on  the  size  of  the  impulse  responses  given  by  Eqn.  (4-1)  is 
not  valid.  As  such  the  impulse  responses  need  not  be  space 
limited  and  may  extend  up  to  the  edge  of  the  array  repre¬ 
senting  the  impulse  responses.  However  as  the  number  of 
transfer  functions  being  multiplexed  increases,  the  size  of 
the  central  array  passed  by  the  low  pass  filter  becomes 
smaller  and  results  in  the  loss  of  more  and  more  components. 
In  general  if  N  transfer  functions  are  being  multiplexed  in 
each  dimension,  the  size  of  the  central  square  array  U  in 
each  dimension  passed  by  the  low  pass  filter  is  given  by 


(a)  Low  pass  filtering  of  transfer  functions:  central 
32  x  32  elements  are  selected  from  each  transfer 
function. 


(b)  Repositioning  of  the  selected  elements  in  the 
composite  array. 


Figure  4-25.  Generation  of  composite  transfer  function 
array  using  low  pass  filtering  technique. 
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o  -  §  ,  (4-2) 

where  b  is  the  size  of  the  composite  transfer  function 
array.  Also,  since  the  elements  from  the  central  portion 
of  the  transfer  function  array  are  repositioned  throughout 
the  composite  array  as  shown  in  Fig.  (4-25)  the  spatial 
frequency  scale  in  the  composite  transfer  function  is  changed 
by  a  factor  of  N  relative  to  the  spatial  frequency  scale  of 
the  individual  transfer  functions.  This  results  in  a  re¬ 
duction  in  the  size  of  the  impulse  responses  during  playback 
by  the  same  factor.  This  multiplexing  scheme  is  suitable  for 
computer  multiplexing  only  as  it  requires  repositioning  of 
elements  and  it  is  not  possible  to  device  a  simple  optical 
equivalent  of  this  multiplexing  technique  for  recording 
the  composite  hologram.  An  experiment  to  verify  this  multi¬ 
plexing  method  was  carried  out  using  the  impulse  responses 
shown  in  Fig.  (4-26) .  Note  that  the  impulse  responses  ex¬ 
tend  upto  the  edges  of  the  array  in  both  x  and  y  directions. 
The  composite  hologram  generated  by  the  computer  is  shown 
in  Fig,  (4-27).  The  result  of  a  computer  simulation  of  the 
playback  of  this  composite  transfer  function  is  shown  in 
Fig.  (4-28) .  Note  that  the  edges  of  the  function  are  not 
sharp  due  to  the  loss  of  high  frequency  components  during 
the  multiplexing  step.  Note  also  that  the  size  of  the  com¬ 
bined  impulse  response  during  playback  is  only  32  x  32  ele¬ 
ments  compared  to  the  original  size  of  64  x  64  elements  as 
a  result  of  scaling  in  the  Fourier  plane. 


multiple 


Figure  4-27,  Composite  hologram  generated  using  low 
pass  filtering  technique. 
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Figure  4-28.  Computer  simulated  output  when  all  the 

impulse  responses  of  Fig.  4-26  are  played 
back. 
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CHAPTER  5 


CONCLUSIONS 

The  representation  of  space-variant  systems  using  en¬ 
coded  reference  beams  requires  diffusers  with  good  correla¬ 
tion  properties.  One  of  the  objectives  of  this  report  has 
been  to  evaluate  the  correlation  properties  of  a  family  of 
binary  phase  codes  for  use  as  diffusers  in  multiplex  holo¬ 
graphy.  In  Chapter  2  the  results  of  extensive  computer  simu¬ 
lations  to  compute  the  autocorrelation  and  the  crosscorrela¬ 
tions  of  a  set  of  codes  described  by  Gold  for  use  in  spread 
spectrum  communication  systems  was  presented.  Simulations 
of  multiplex  holography  using  Gold  codes  of  different  lengths 
were  also  carried  out.  The  results  of  these  simulations  indi¬ 
cate  that,  due  to  non  ideal  correlation  properties  of  the  Gold 
codes,  the  magnitudes  of  the  crosstalk  terms  are  quite  large 
resulting  in  poor  reconstruction.  However  it  was  observed 
that  when  the  impulse  responses  have  very  small  spatial  widths, 
acceptable  levels  of  the  signal-to-crosstalk  ratio  were  ob¬ 
tained.  Thus  it  may  be  concluded  that  the  method  of  space- 
variant  system  representation  using  the  Gold  codes  as  dif¬ 
fusers  is  more  suitable  for  applications  involving  space- 
variant  systems  such  as  magnifiers  which  transform  points  in 
the  input  plane  to  points  in  the  output  plane,  resulting  in 
delta  function  like  impulse  responses. 

The  effect  of  chirped  wave  illumination  was  briefly  des¬ 
cribed  and  the  need  for  fabricating  near  perfect  phase  dif- 
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fusers  was  also  demonstrated.  A  technique  for  fabricating 
phase  masks  using  dichromated  gelatin  is  discussed  in  Appen¬ 
dix  G.  Additional  work  needs  to  be  done  in  this  area  to  perfect 
this  fabrication  process  as  well  as  to  determine  simpler  meas¬ 
uring  techniques  for  evaluating  the  quality  of  the  resultant 
phase  masks. 

Another  objective  of  this  research  was  to  develop  an  al¬ 
ternative  multiplexing  technique  to  generate  composite  holo¬ 
grams  representing  the  system  transfer  functions.  In  Chapter 
3  a  technique  in  which  the  transfer  functions  are  sampled  in 
the  Fourier  plane  and  repositioned  to  represent  a  composite 
hologram  was  presented.  The  multiplexed  hologram  generated 
by  this  technique  contains  samples  of  the  transfer  functions 
in  nonoverlapping  regions  and  hence  the  problem  of  hologram- 
to-hologram  crosstalk  is  completly  eliminated.  The  method 
also  requires  a  single  reference  beam,  unlike  the  encoded  ref¬ 
erence  beam  approach  that  required  a  number  of  reference  beams. 
However  this  technique  requires  generation  of  multiple  images 
of  the  input  function  during  the  playback  step.  Several  schemes 
which  permit  generation  of  these  coherent  multiple  images  were 
briefly  described.  Additional  research  to  implement  these  or 
other  methods  for  multiple  imaging  needs  to  be  done  in  order 
for  this  multiplexing  technique  to  become  practicable  for  the 
representation  of  space-variant  systems  characterized  by  a 
large  number  of  impulse  responses. 
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Experimental  results  using  computer  multiplexed  holograms 
to  represent  a  space-variant  system  sampled  at  2  x  2  points 
in  the  input  plane  were  presented  in  Chapter  4.  The  experi¬ 
ments  were  conducted  for  both  disjoint  impulse  responses  as 
well  as  for  overlapping  impulse  responses.  The  property  of 
coherent  addition  that  is  required  in  the  case  of  overlapping 
impulse  responses  was  also  verified  through  these  experiments. 

A  slight  variation  of  this  multiplexing  technique  using  a  low 
pass  filter  in  the  transfer  function  plane  followed  by  repo¬ 
sitioning  of  the  filtered  components  was  also  presented.  A 
combination  of  these  techniques  may  be  adopted  to  multiplex 
larger  number  of  transfer  functions  in  a  single  composite 
array.  Multiplication  of  the  impulse  responses  by  a  random 
phase  mask  to  distribute  the  transfer  functions  more  evenly 
so  as  to  reduce  the  quantization  losses  of  small  components 
during  the  generation  of  computer  multiplexed  holograms  was 
also  demonstrated.  The  computer  simulations  and  the  experi¬ 
mental  results  presented  in  this  report  demonstrate  the 
ability  of  the  sampled  input/sampled  transfer  function  approach 
to  effectively  represent  any  slowly  varying,  linear,  space- 
variant  system  with  finite  spatial  extent  of  the  impulse  re¬ 
sponses.  However  implementation  of  this  method  for  very  large 
size  sampling  arrays  in  the  input  plane  requires  high  precis¬ 
ion  in  the  alignment  of  the  multiple  images  of  the  input  func¬ 
tion  during  playback  and  is  likely  to  be  a  limiting  factor 
in  practical  systems. 
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APPENDIX  A 

Computer  Program  to  Generate  Gold  Codes  of  Length 
511  Bits  and  a  Set  of  Nine  Codes  Generated  by  the  Program. 
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PROGRAM  COOt 

This  PROGRAM  GENERATES  a  set  of  the  golu  code 
SEQUENCES  OF  LENGTH  512  BITS  bACH.  THE  REQUIRED 
I  N  0  U  T  S  TO  THE  PROGRAM  ARE  (1)  THE  SHIFT  REGISTER 
CONNECTIONS  REPRESENTING  THE  POLYNOMIAL/  (2)  THE 
'UMBER  OF  CODES  IN  THE  SET  (LIMITED  TO  A  MAXIMUM 
UF  25.)  THE  PROGRAM  USES  000000000000000001  AS 

the  first  seed  for  generating  the  first  code  in 
The  set.  SUBSEQUENTLY  the  VaLU>=  of  the  seed  is 
incremented  ry  i  and  the  new  seed  is  checked 
tu  verify  whether  it  is  a  segmfnt  of  the  codes 

ALREADY  GENERATzD.  if  SO  THE  SEED  IS  AGAIN 

incremented  py  i  and  again  checked,  the  output 

of  THE  PROGRAM  IS  A  SET  OF  UNIQUE  CODES.  THE 
RROGRAM  PRINTS  THE  SFED  USED  AND  THE  corresponding 
HOLD  CODE  GENERATED. 

DIMENSION  I Z ( 25 / 530 ) / I  A  C 1 8  > 

5  FORMAT<////5 

READ  THE  SHIFT  REGISTER  CONNECTIONS, 
bH IFT  REGISTER  CONNECTIONS  ARE  ENTERED  ON  THE 
datacard  starting  from  the  highest  order. 

For  EXAMPLE  1+X+X**3*X**5+X**6+X**8+X**11+X**12+X**15 
*X**16+X**18  IS  ENTERED  AS  1 01 1  0  011  i)  01 01101 0 1 
IGNORING  THE  CONSTANT  1, 

REAJ(5,10) ( IA( I), 1=1/18) 

10  FORMAT(IBU) 

••RITE  (6/ 5) 

•  nITE<6,15){Ia<I>, 1*1,18) 

15  F  0RMAT(5X,26HSHlFT  REGISTER  CONNECT  I  OinSz /5X,  16  1 1, //// ) 
HEAD  THE  DESIRED  NUMBER  OF  CODES  IV  THE  SET, 
■<EAD(5/20)N 
20  F o R M  A  T  < I  3 ) 

•>RITE(6,25)N 

25  FoPHAT(5X, 20HUU1PER  OF  SEQUENCES: , I  3, //// ) 

C  SET  THE  SEED  FOR  THE  FIRST  CODE  AS  0  0  0  0  U  0  0  0  0  0  Ij  0  0  0  0  0  0 1  , 

i'03i)  J*l,  17 
30  IZ(X/J)*0 
IZ<1,18>*1 
■'0300  K s  1  /  N 
lF(K.EO,l)P0,5f 

C  iNCRcME'N’-  THE  '  -D  USED  FOR  THE  =>REVIUUS  CODE  «Y  1. 

50  liQ52  J  =  i  3 
52  IZ(K/J)=IZ(K-1,J) 

55  DO  6 0 Js 1 , 18 

I F ( ( I Z ( K , 1 9- J ) ) ,FO.l)60,56 

56  IZ(K,19-J>=1 
I J= J-l 

l' 0 5 8  JJ=1,IJ 
58  I Z ( K / 19- JJ ) sO 
G0T062 
60  CONTINUE 
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VERIFY  WHETHER  THE  NEW  SEED  IS  A  SEGMENT  OF  THF 
3REVI0USLY  GENERATED  CODES  IN  THE  SET ,  IF  SO  RETURN 
TO  THE  PREVIOUS  STEP  AND  CHANGE  THE  VALUE  OF  THE 
SEED.  IF  NOT  USE  THIS  SEED  FOR  GENERATING 
ANOTHER  WEM8ER  IN  THE  SET, 

62  00701*1,512 
KKsK-1 
U07UI*1,KK 
D065J*1, 18 

IF( ( IZ(K»19-J> ) ,EQ. ( IZ( I , 18+L- J) ) >65,70 
65  CONTINUE 
G  0  T  0  5  5 
70  CONTINUE 

80  R  I T  E  (  6  ,  3  5  )  (  IZ(^,J),J*1,16) 

85  FqRHAT<5X, 21HSEED  FQ9  THE  SEQUENCE , /5X , 18 1 1 ,//// > 

C  GENERATE  THE  GOLD  CODE  USING  THE  SEED 

C  PREVIOUSLY  SELECTED, 

DO150 I =19 , 530 
“1  =  0 

L-U90  J=l,  18 

90  f =Im( J)*IZ(K, J*I-19)*M 
RsH/2. 0 
N  M  a  R 
W=R-MM 

150  1Z(K,  I  )=R*2. 0 
*RITS(6,1S0> 

180  r ORH4T<5X,26HGOwD  CODE  FOR  THIS  SEED  IS,//) 
‘■RITS(6,2uO)  ( IZ(K,  J),  J=lo,530) 

200  F0RMAT(5x,48U  ) 

3  U  0  svRITS(S,3) 
end 


! 
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APPENDIX  B 

Computer  Program  to  Evaluate  the 
Correlation  Properties  of  the  Gold  Codes 
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PROGRAM  SP ACEV  A  H 

this  program  computes  and  plots  the  autocorrelation 

DP  A  GOLD  CODE  IN'  A  SET  OF  NINE  128  BIT  CODES  AND 
ITS  CROSSCOHRELaTION  WITH  THE  OTHER  EIGHT  MEMBERS 
IN  the  SET.  THIS  IS  DONE  by  MULTIPLYING  THE  FOURIER 
transforms  of  the  codes  and  then  fourier  transforming 
the  PRODUCT.  IN  THIS  PROGRAM  ONLY  the  CENTRAL  42 
SITS  OF  THE  128  BIT  CODE  ARE  USED  IN  THE  CALCULATION 
SO  THAT  THIS  OUTPUT  KAY  BE  COMPARED  WITH  THE  OUTPUT 
OF  THE  PROGRAM  HPXHOLO. 

DIMENSION  A(l28).G(12U)»H(128)fT(l28)  ,B(9,l28) 
ty°e  complex  G. h.t.cmplx.conjg 
callplotsi  n,  (i,i) 

WRITE! 6, 220) 
no 5 1  si, 12 ti 
H( I)S(0, 0,0,0) 

5  T( I )S( 0.0, n.0) 

READ  the  12a  81 T  CODE  TO  BE  USED  AS  THE  COMMON  CODE 
In  THE  CALCULATION  OF  CORRELATION  WITH  A  SET  OF  NINE 
129  BIT  CODES,  the  code  is  read  In  as  ZEROS  AND  ThOS, 
RE AD( 5, 10 ) { A( I ) , I s 1, 128 ) 
io  form at< 64F1 , 0) 

CONVERT  THE  VALUES  IN  THE  CODE  TO  *1  AND  -1  TO 
REPRESENT  A  PHASE  MASK  WITH  180  DEGREES  PHASE 

difference. 

•’0201  =  44,  85 
I F ( A ( I ) , GT . 0 ) 15, 12 
12  »(I)*-lt0 
30  TO  ?0 
15  M I )=1 
20  CONTINUE 
00301=44,65 
“( I )=CMPLX(A<  I), 0,0) 

30  CONTINUE 

calculate  the  Fourier  transfohm  of  the  code  . 

CALL  FOURIER  <H> 

calculate  the  conjugate  of  the  fourier  transform  of 
the  code  and  stjre  in  the  array  t, 

JO  40  1=1,128 
40  T ( I ) =C0N JO ( H ( I ) ) 

The  FOLLOWING  dj  LOOP  READS  all  the  nine  codes  IN 
the  set  one  at  a  time  and  computes  their  correlation 
<■  I  Trt  the  code  prfviusly  read,  this  set  of  NINE  CODES 
*LS0  INCLUDES  The  PREVIOUSLY  READ  CODE  as  A  member 
ANO  HENCE  ONE  of  THE  OUTPUTS  IS  the  AUTOCORRELATION 
and  the  REST  AR=  THE  CPOSSCORRELAT I ONS , 

DG12Q  KK s 1 , 9 
DO  60  1=1,128 
G(I)=(0, 0,0,0) 

60  h ( I)=(0, 0,0.0) 


no  ci  on  onoo  oc> 
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*E A L)  A  MEMBER  FROM  THE  SET  AWu  CONVeRT  THE  VALUES 
TO  *1  A,\P  -1, 

HEAU<5,10) <A< I ) , 1=1,128) 

61  r0RMAT(5X,64Fl.U) 

»R I  TEC  6*  61 ) ( A ( I ) , 1=1,127) 
vRITE(6,220) 

DO  70  1=44,85 
IF(A( !) ,GT. 0)65,62 

62  A ( i ) s • 1 
GO  TO  70 

65  A ( I ) s 1 
70  CONTINUE 

l>0  80  1=44,85 
M(I)=CMPLX(A(I),0,0) 

80  CONTINUE 

calculate  the  fjukjer  transform, 

CALLFOURIER(H) 

MULTIPLY  THE  FOjPISR  TRANSFORM  OF  THE  CODES, 

DO  V0  1=1,128 
V0  3( I )=T< I )*H( I ) 

COMPUTE  THE  FOURIER  TRANSFORM  OF  THE  PRODUCT  TO  OuTaIn 
AS  THE  OUTPUT  THE  CORRELATION  BETWEEN  THE  CODES. 

call  Fourier  <g> 

DO  120  1=1,128 
120  F  t  4K , 1 )  =  CABS ( G(  I)) 

FORMALIZE  the  MAGNITUDES  of  all  the  elements  in  the 
UUTPJTS  WITH  REFERENCE  TO  THE  VALUE  OF  THE  LARGEST 
ELEMENT  IN  THE  =NTIRE  SET, 

XM*0,U 
'0122  1=1,9 
u0122J=l,123 
122  XM=MAX1F(B( I, J),XM) 

001251=1,9 
50125J  =  1, 128 
125  H( I, J)=B( I , J ) *9 V/X M 

plot  each  of  th=  normalized  outputs  to  a  width  of 
'8,56  INCHES, 
uO?  0  OK  =  1 , 3 
U0400L*1,3 

CaLLpLOT(0,0,-2,56,2) 

C ALLpL0T ( 0 , 0, 0 , J, 2) 

J=(K-1)*3+L 
L:  01601=1, 128 

160  CALLpL0T(b( J, I)/40,-I*0.02,2) 

C  PRINT  THE  NORMALIZED  VALUES  OF  THE  CORRELATION  OUTPUT, 

.RITE(6,2U0)(B(J, I), 1=1,128) 

200  FORy,AT(5X,  16F4,  j) 
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'RITE<8,220) 

220  FOR.IAK//) 

400  CALLPL9T(0.U,-3.0,-3) 
500  CaLLPL0T<3.G,9, 0,-3) 
CALL°LOT(0,0, 0, 0,999) 


^URROUTINfc  FOURIER  («) 

this  subroutine  calculates  the  foupier  transform 

(F  AN  ARRAY  CF  126  ELEMENTS, 

this  program  also  shifts  the  elements  in  the 
array  to  take  care  of  the  fft  algorithm  which 
ASSUMES  the  first  element  as  the  origin, 
dimension  M(3),S(32), INV(32)  ,R<128) 
type  cokplexx.b 

JATA(Mr7, 0, 0) 

I  T  =  U 

10  io  20  1=1,84 
■A  =  3  (  I  ) 

*<!)•  5(1+84) 

•’(  I  +64  )  =X 
20  CONTINUE 

I F ( IT. EO. 1)30, 22 

c  call  fft  algorithm  harm, 

22  call  harm  <8,m, inv.s, i, IFErr) 

IT  =  1 

iO  TO  1U 
30  RETURN 

end 


APPENDIX  C 

Computer  Program  to  Simulate  the  Output 
of  a  Multiplexed  Hologram  Using  Gold  Codes 
as  Phase  Diffusers  in  the  Reference  Beam  Path. 
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PROGRAM  MPXHOLn 

C  This  PRuGKAM  SIMULATES  Th£  output  of  A  SYSTEM 
C  -SING  The  GOLD  CODES  AS  D I FFUSER  5  I  N  Mb 

C  -<EFERENCE  REAM  3 ATH  TO  REPRESENT  A  SPACc-V  AR I  AnT 

C  PROCESSOR. A  TOTAL  OF  TWO  TRANSFER  FUNCTIONS 

0  '‘E=»RESENTIK'G  TH  =  SPACE-VARIaNT  SYSTEM  ARE  MULTIPLtXED 

C  In  A  SINGLE  ARRAY.  THE  PROGRAM  ALSO  PLOTS  THE 

C  I  m»uLSE  RESPOnSzS  USED  IN  THE  SIMULATION  AND  TRt 

C  OUTPUT  WHEN  THEY  ARE  RECORDED  AND  PLAYEu  HAC* 

C  ‘jSInG  Trie  GOLD  CODES  WITHOUT  ANY  MULTIPLEXING, 

c  in  this  program  only  the  central  42  pits  of  the 

C  12*  R I T  COOP  A R r  USED  IN  THE  CALCULATION  SO  TH»T  THt 

C  UTPUT  OF  THE  SIMULATION  WITH  THE  TOTAL  NINTH 

c  equal  to  the  sum  uf  the  widths  of  the  t*p  cupes 

C  AND  THE  WIDTH  0**  T Hg  IMPULSE  RESPONSE  IS  LESS  Tm A  j 

C  The  SIZE  of  the  OUTPUT  ARRAY,  FOR  THE  SAME 

C  REASON  THE  MPU.SE  RES3OVSES  ARE  ALsO  L  M  i  TeO 

C  To  A  SIZE  OF  42  HITS. 

DIMENSION  A(128),G<128),H<128).S(l2d),R(12'>),T(i28) 
Type  COMPLEX  3,m,T,CMPLX,C0NJ3  ,s 
0ALLdL0T3< 0, o,  1 ) 

-'0  2  1=1, 12a 

5  1  <  I  >  S<  n.  o,  II.  0  ) 

*-ulOOK  =  l,2 
'■’071=1,128 
7  3( n  =  u.o 

^E AD  A  12a  BIT  GOLD  CODE,  VALUES  READ  IN 
«k=  ZEROS  AND  twos. 
h£AD(5,10) (  A  (  I  )  , 1=1,128) 

10  F  UR.1AT(64F1,  0) 

cowert  the  values  to  *i  and  -i  to  represent 
A  =hAS=  Mask  with  isn  DEGREES  PHASE  DIFFERENCE 
"ETwEtN  THE  ELEMENTS. 

JO  20  1=44,65 
I  F  {  A (  I  ) , GT , 0 )  1?,12 
12  M  I  )*«1 
*•0  TO  20 
15  A ( I  )  =  i 
20  CunTIWUE 
'■'0  2p  1  =  1,128 
H(  i )  =  < o , n, n , o ) 

25  H  (  l)  =  (0, 0,0,0) 
i'O  GO  1  =  44,85 
■*  <  I  )=CMPLX<  a{  I  )  ,  0  .  0  ) 

30  CONTINUE 

Fourier  transform  the  array  to  represent 

The  REFERENCE  REAM  ILLUMINATING  THE  HOLOGRAM, 

Gall  FOURIER  (H) 

00321=1,128 
32  A (  I  )  s o  ,  o 


o  o  o  o  o  o  o  o  n  n  o  o  on 
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C  read  a m  I ;,DuiSt  response  yf  The  ^pace-vaR  I  a,»t  systpim. 

'EAD(5,33)  (  A  (  J  )  ,  Js44»8=5) 

33  F  CR.MAT(21F2,0) 

:l035  1*1,128 
‘(!)=C*1PLX(A(I),0.0) 

35  R<  I ) =A?S ( A (  I  )  ) 

c  plot  and  print  the  impulse  responses. 

CALLPLOTT(B) 

CaLL®L0T(3.0, 0. J,-3) 

C  COMPUTE  THE  TRANSFER  FUNCTION  , 

CALL  FOURIER  (G) 

•’E'IERaTE  THE  COMPOSITE  TRAMSFfeR  ^UNCTION  ARRAY  tJY 
SUMMING  THE  PROJUCT  OF  THE  FOuRIER  TRANSFORMS 
OF  EACH  OF  THE  IMPULSE  RESPONSES  AND  THE 
CORRESPONDING  REFERENCE  BEAM  FUNCTION, 

OQ  40  1=1,128 

40  T(  I )sTC I WG( I )*J0NJG(H( I ) ) 

SIMULATE  THE  PLAVbAC*  OF  AN  IMPULSE  RESPONSE 
»H=N  THE  GOLD  CODE  IS  USED  IN  THE  RECORDING 
ANO  THE  PLAYBACK  STEP. 

D05ul=l,12fl 

50  U< l )sH< I )*CONJG(H(I ) )*S( I  ) 

call  Fourier  (G> 

nO  60  1=1,128 
60  -< I ) =C AbS ( G (  I  )  ) 

PLOT  AND  PRINT  |Hb  OUTPUT  WHEN  THE  IMPULSE  RESPONSES 
ARE  RECORDED  ANj  PLAYED  RACK  USivG  THE  COLD  CODES. 
CALL3LOTT(?) 

CALLBLOT(-3. 0,-4, 0,-3) 
lF(K.tH.l)70,10 J 
70  f>0*U  1=1,126 
60  S ( I ) s H (  I  ) 

100  CONTINUE 

CALL3L0T(6, 0,8. 0,-3) 

-'01131*1, l?b 

SIMULATE  THE  PLaYOACK  OF  THE  SYSTEM  WH  =  h  THE 
TRANSFER  FUNCTIONS  ARE  ACCESSED  (A)  INDIVIDUALLY 
anD  (rf)  SIMULTANEOUSLY  FROM  ThE  COMPOSITE  ARRAY. 

■’<  I  )=T(  I  >*S<  I  ) 

T  (  1 ) =T  < I ) *H  C  I  ) 

110  s ( I ) =  G ( I ) ♦ T  <  I  ) 

CALL" UURIER(G) 

CALLFOURlcR(S) 

CaLLFOJRIERI  T  ) 

•)U  1201=1, 128 
120  n(I)*CAbS(C( 1 )) 

PRI.mT  and  PLOT  THE  OUTPUT  WHEN  THE  FIRST  TRANSFER 
FUNCTION  IS  ACCESSED  FROM  THE  COMPOSITE  ARRAY, 
CALL®LOTT(:R) 

CALLpL0T(0. 0,-4, 0,-3) 
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.'ul3fiI=l,x2o 
140  !*(  I  )sCAuS(T(  I  )  ) 

C  3 P!  I  i\ T  AND  PLOT  1  He  OUT°UT  N  H  E N  THE  SECOND  T HAnSFER 

C  FUNCTION  IS  ACCzSSED  FROM  THE  COMPOSITE  ARRAY. 

CALL°LOTT(R) 

CaLL°LOT(3.0,2.U,-3) 

001401 =1,126 
140  H(I)sCAUS(S(I>) 

Z  PRInT  and  PLOT  THE  OUTPUT  NHEn  SOTh  THE  TPmvSPER 

C  Fu  MOTIONS  ARE  ACCESSED  FPCM  ThE  CO-'POSITE  ARRaY. 

cALL3LOTT ( H ) 

CALLPLOT(0, a.99»> 
h  n  0 


SURROuTlNt  FOURIER  ( R  > 

c  this  subroutine  calculates  the  fdurifr  t=a  'sform 

C  1»F  AM  array  OF  128  ELEMENTS. 

C  This  PROGRAM  aLbO  SHIFTS  the  ELEMENTS  In  the 
C  “SPAY  TU  TAKE  C*PE  OF  THE  FFT  ALGORITHM  UhIQH 

z  s s u ^ e s  the  first  element  as  the  origin, 

•  ’  I m  E  M  3  I u  M  M ( 3 ) , 3 ( 3  2 ) # I  N  V  (  3  2 )  ,3(128) 

T  Y  3E  UOUPLtXX  ,  R 
ATA(U  =  7,  (i,  0) 

I  T  =  0 

10  i0  20  1=1, M 
>*P( I ) 

3(1)=  0(1+6**) 

--(  1  +  64  )  =  X 
20  00  MT  I  N'JE 

IF( IT, SO. 1)30, 22 
-  -ALL  FFT  ALGORITHM  HARM, 

22  Gall  harm  (o,M, IMV.S.l, IFERR) 

I  T  =  1 

JO  TO  1U 

30  RETURN 

end 
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''U,Rf)UTI\b  PI.UT7CB) 

this  s'Jbpuuti.mf  prints  and  pluts  the  maunj tude*5 

i’F  THE  OUTPUTS,  THE  OUTPUTS  ARE  NORMALIZED  4 1 T  H 
REFERENCE  TO  THE  LARGEST  ELEM£MT  IN  THE  ARRAY. 
DIMENSION  B { 1 2b ) 

<M  =  U  ,  0 
•‘OlUJsl  ,128 
10  <K3MAXlF(d(J)»XH) 

^OPOJsl, 128 
20  »'(  J)sd(  J)*99/XM 

CaLLPLOT(U.O,-?.?6,2) 

CallPL'1T(U.0,0.  J,2) 
ii030I«4*52B 

SO  -ALL3LOT(bM  )/4j,-I*0.  P2,  2) 

-  RITEI6, 40) (E( I > , 1=1, l2b> 

40  FGRmAT<5»X,  32F3.  J) 

*R I TE ( 6, 50) 
pO  FORMAT!////) 

CALL3LOT<C.U,0. J »  3 ) 

RETURN 
Eh  D 


APPENDIX  D 

Computer  Program  to  Evaluate  the  Correlation  Properties 
of  the  Gold  Codes  Illuminated  by  a  Spherical  Wavefront 
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C  SPHERICAL  WAVEFRONT.  THE  COMPUTATIONS  A HE  CARRIED 

C  OUT  FOP  DIFFERENT  VALUES  OF  THE  CHIRP  AS  SPECIFIED 

C  eY  THE  RADIUS  OF  CURVATURE  OF  THE  WAVEFRONT  AND  THE 

C  r:  I OTH  OF  THE  mask,  THE  OUTPUTS  ARE  FORMALIZED  r  Y 

C  FORCING  THE  AREA  UNDER  EACH  OF  THE  AUTOCORRELATION 

C  PEAK  FOR  DIFFERENT  VALUES  OF  THE  CHIRP  TO  RE  EQUAL 

C  SO  THAT  THE  OUTPUTS  MAY  PE  COmPa3ED  WITH  EACH  OThEH. 

C  AN  ESTIMATE  OF  THE  NOISE  TO  SIGNAL  RATIO  IS  ALSU 

C  ''ADE  d Y  CALCULATING  THE  RATIO  OF  A»EA  UNDER  THE 

C  CORRELATION  CUR  YES  TO  THE  AREA  OF  THE  AUTOCORRELATION 

C  PEAK. 

DIMENSION  At128),H(2D49),T(204S>,ytl28).Ctl29>.D(l2«> 
T  Y3E  COMPLEX  H, T,CMPLX,CONJG,C,U, ARA,.  RC 
CALLRLOTStfl, 0, i> 

C  READ  THE  TWO  128  HIT  CODES.  CONVERT  THE  VALUES  TO 

C  +1  AND  -1  TO  REPRESENT  A  PHASE  MaSK  WITH  180 

C  JtGREES  phase  difference. 

REAU(5,10) { A< I ) , I  si, 128 >  •(&<!)• 1 >1# 120) 

in  FORHATt 64F1, 0) 

00  2n  1=44,68 
I  Ft  At  I )  ,  OT , 0 )  13,12 

12  A(  I)=-l 

GO  TO  20 
15  A ( I ) s 1 
20  CONTINUE 
l-'0?6l  =  44,  b5 
lF(Ut I). GT, 0.0)25, 22 
22  "  (  I  )  =  - 1 
-0T026 

25  n(  I  ) =i 

26  CONTINUE 

C  IhE  FOLLOWING  DO  LOOP  COMPUTES  T H£  AUTOOORPELAT I  ON 

.C  OF  (A)  WITH  ITSrLF  AND  CROSSCURREL  AT  I  ON  OF  (A)  WITH 

C  <B)  FOR  DIFFERENT  VALUES  OF  CHIRP, 

L‘O5u0KK*l,  8 
0051=1,2048 

h  ( I )  =  ( o .  o ,  n .  o  > 

5  Tt I )  =  <  0, 0, 0,  0) 

C  "EAC  THE  RADIUS  OF  CURVATURE  OF  THE  SPHtRlCAL 

C  •  AVEFROnT  AND  THE  WIDTH  OF  THE  CODE  MASK. 

''EADt5,d)R,W 
9  FoRMATf 2F6.2) 

C  - 0 A L)  THE  CODES  IN  A  LARGER  ARRAY  SO  THaT  EACH  ELEMENT 

C  IN  THE  ORIGINAL  ARRAY  OCCUPIES  16  ELEMENTS  IN  THE 

C  NEW  ARRAY,  THIS  IS  NECESSARY  TO  REPRESENT  THE  PHASE 

C  VARIATIONS  WITHIN  EACH  ELEMENT  WHEN  THE  CODE 

C  IS  ILLUMINATED  8Y  A  SPHERICAL  WAVEFRONT, 


no  u  u  o  i  j  u  u 
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1-0  30  1  =  4 *,&5 
l‘O30J*1#16 

T( ( I-1)*16*J)=CkPLX(«< I ) ,0.0) 

30  M(  ( I-1)*16*J)=CmPLX(  A  (  I  >  #  0  •  0*> 

C  CALCULATE  THE  PATH  DIFFERENCE  AT  THE  CENTER  OF  EACH 

C  cL?HENT  RELATIVE  TO  THE  CENTER  OF  THF  ARRAY, 

P03*>  1=1,1024 

S= ( SORT  <R**2*(<I-0,5)*w/2049)**?)-R) 

COMPUTE  THE  PHASE  DIFFERENCE  FOR  A  WAVELENGTH  OF 
300  NANOMETERS, 
i-'sS*(in**4)/5 
L  =  D 

i's(D-L)*2*3,142s 

change  the  values  of  the  elements  in  the  complex  array 

To  INCLUDE  THE  ;FFECT  OF  THE  PHASE  CHANGE  DUE  TO 

the  curvature  of  the  wavffront. 

T (1024*1 )=T (1024*1 )*CMPLX<COS(D> ,SIN(D) ) 

T( 1024-1 )=T( 1024- I )«CMPLX( COS( D) ,SIn(U) ) 

M( 1024*1 >*H(J 024*1 >*CMPLX<COS(D> #SIN(D) ) 

35  h (  1024-1  )=H( 102f- I >*CMPLX< COS  CD), SIN(D) ) 

COMPUTE  THE  AUTj  AND  THE  CROSS-CORRELATIONS, 

CALL  FOURIER  <H) 

CALLrOURItR(T) 

CO  40  1=1,2046 
T ( I ) =  T ( I )*COMJG(H(  I)) 

40  *<  I  )=H(  I  )  *CONi  JG  (  H  (  I  )  ) 

CALLFOUHlcR(T) 

call  Fourier  <h) 

COMPUTE  THE  AVERAGE  MAGNITUDE  OF  EVERY  16  ELEMENTS 
4nT  GENERATE  OUTPUT  ARRAYS  OF  12?  ELEMENTS  EACH, 

00  120  1=1,129 
U=< I-l)*lft*l 
M  =  J  *  1 4 
DO110K=J,M 

T ( J)=T( J)*T(K*1) 

110  *( J)=H( J)*H(K*1) 

D  (  I )  =  T  (  J ) 

120  C( I )=H( J) 

C  COMPUTE  THE  MAGNITUDE  OF  THE  PEAK  OF  AUTO-CORRELATION, 

X  M  =  0 , 0 

130  '>01.351  =  1, 126 
X=CARS(C<  I  )  ) 

135  X,M  =  MAX1F  (X,XM 

C  COMPUTE  THE  ALGEBRA  I C  SUM  OF  THE  ELEMENTS  IN  THE 

C  OUTPUT  ARRAYS. 

ARA=<0, 0,0.0) 

ARC=(0, 0,0,0) 

140  D0150  1=1,120 

ara  =  ara  +  C(  I  ) 

150  ARC  =  ARC*D(  I  ) 

AA=CABS( ARA) 


i k 


nooon  (H)uo(i  ooo 


*CsCAbS( AkC) 

COMPJTE  ThF  NOISE  TO  SIGNAL  RATIO  OF  AUTOCORRELATION 
AND  THE  CKOSSCORRELATION  FUNCTIONS  USING  THE  PEAK 
CF  AUTOCORRELATION  AS  REFERENCE. 

WA= A A/XM-1 
RC=AC/XM 

formalize  the  values  with  reference  to  the  peak  of 

-UT0C0RRELATI0N. ALSO  SCALE  THE  VALUES  ACCORDING  TO  Tnfc 
MOTH  OF  THE  MASKS  WITH  WlOTHs3.Bl  AS  REFERENCE, 

this  ensures  that  the  area  under  the  autocorrelation 

PEAKS  FOR  MASKS  with  DIFFERENT  WIDTHS  ARE  ALL  EQUAL. 
150170  1*1,129 

D(  I  )*D<  I  )*{99/X'1>*<3 ,81/W) 

170  C< I  )=C( I )*<99/XM  W3.81/W) 

• RITE(6,190)R 

190  F0RiAT(5X,21HRA JIUS  OF  WAVEFRONT  =,F6,2) 

WRITE(6,200)W 

200  F0RMAT(5X,21HwlDTH  OF  CODE  MASK  *  ,F6,2) 

■••■RITE  <6, 210  )XM 

210  ► 0RMAT(5X,31WHEIGHT  OF  AUTOCORRELATION  PE AK = , E9 , 2// > 
w RITE(6,220)AA 

220  F0Rr',AT<Sx,2lHAAsAREA  OF  auTOCqRN  ,E9,2) 

■RITE<6,230) AC 

230  FOR 1 1  AT (5x,21uACsAR5A  OF  CRDSSCORN  ,E9,2) 

..RITS(6,240  )RA 
240  *■  O R n A T  (  5 X  , 

14HNDISE  TO  5  I G  \AL  RATIO  OF  A JTOCORRELAT I  ON* , F8 . 2//> 
'*RITE(6,250)RC 
250  •" 0 R i . A T  ( 5 X , 

15-1HRATI0  OF  NOISE  OF  CROSS  CORN  TO  SIGNAL  OF  AUTOCORNs# 
2FS.2//) 

plot  the  auto  a \jd  the  cross-correlations.  the  width 

OF  PLOTS  APE  SCALED  ACCORDING  TO  THE  WIDTH  OF  EACH 
'VASK.  <  WIDTH  OF  aLUTS  -FOR  MASK  WITH  W=3.61  IS 
TAKEN  as  2,56  IviCHES  A''D  is  USED  AS  reference  FOR 
COMPUTING  THE  P..OT  SIZE  FOR  OTHER  MASK  wlDTHS.) 
CALL:,L0T(U.C,-2,56*W/3,81,2) 

CALLoL0T(0.0,0.i),2) 

003001*1,128 
x  =  C43S(C(  I  )  ) 

30  0  C ALLPLOT ( X/100, -1*0, 02* w/3,dl* 2) 

CALLPLOKO,  0,-4, 0,-3) 

CALLaLOT(0.0,-2.56*W/3,Hl,2) 

:allrlot(o.o,o. 0,2) 

004001=1,126 
XsCA3S(D( I ) ) 

400  CALLPLOT(X/ino, -1*0,02* W/3, 81,2) 

CALL3L0T(3, 00, 4, 00,-3) 

500  CONTINUE 

CALLRL0T( 0 . 0, 0,0,999) 

End 


SlHKOUTl'Jfc  FOURIER  (?) 

iJ  I  M£MS  1 1)\'  M(?),  5(512),  I  n|V(5l2)*B(  2048) 
T Y3E  COi*lPLEXX#R 
!'<ATA(M*11,  0*  0  ) 

I  T*0 

10  DO  20  1*1,1024 
**R<  I  ) 

h(l)s  ■<(  I  +1  £3  24  ) 

*( I+1024)=X 
20  CONTINUE 

IF( IT. £0,1)30, 22 

22  "ALL  HARM  (B,M, 1  MV, S, 1, I FERR ) 

I  T  *  1 

'<0  TO  10 
30  «ETURN 
cNO 


APPENDIX  E 

Computer  Program  to  Simulate  the  Output 
of  a  1-D  Processor  Using  the  Sampled  Transfer 
Function  Approach  for  Multiplex  Holography 
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PROGRAM  SPCEHIV^ 

This  PROGRAM  MULTIPLEXES  the  TRANSFER  FUNCTIONS 
OF  FOUR  IMPULSE  RESPONSES  AND  GENERATES  A 
SINGLE  COMPOSITE  ARRAY  USING  THE  SAMPLING 
TECHNIQUE  IN  THE  TRANSFER  FUNCTION  PLANc. 

the  program  plots  the  inputs  representing  THE 
impulse  response  and  their  respective 
transfer  functions  before  and  after  sampling, 
the  program  also  simulates  and  plots  the 

PLAY3ACK  OF  THE  IMPULSE  RESPONSES  WHEN 
The  TRANSFER  FUNCTIONS  ARE  ACCESSED  (A)  INDIVIDUALLY, 
<B)  simultanfos_y. 

I  MENS  I  ON  A(l2R),B(l2b),G(128),H(l28),T(128> ,C(12B) 

type  complex  g,m,t,cmplx,conjg,x 
call  plots!  o, o, i> 

POlOOKsl,  4 
C  READ  TME  IMPULSE  RESPONSE. 

READ(5,10HA(l),Isl,12S) 

10  FORMAT(64Fl, 0) 

00201=1, 128 

20  H( I)*C4PLX( A( l  ) ,  0, 0) 

C  GENERATE  THE  TRANSFER  FUNCTION  USING  THE  FAST 

C  fojhier  TRANSFORM  ROUTINE. 

Call  FOURIER  (H) 

c  sample  the  transfer  function  at  an  interval  of  four 
c  elements  and  load  the  SAMPLES  IN  THE  COMPOSITE  array. 
003U  I=K,128,4 
30  TCI)sH(I) 

C  SCALE  THE  IMPULSE  RESPONSES  TO  A  MAXIMUM 

C  VALUE  OF  99  AND  PRINT. 

XM  =  0 , 0 

L’O  40  1=1, 12P 
A(  I ) = A  RS ( A  (  I  )  ) 

40  XM  =  MiAX1F(A(I),XM) 

00501=1,128 
50  A { I ) s  A ( 1 ) *  99/ XM 

'RITE!6,  60) <A< 1), 1=1,128) 

60  F0RMAT(5X,32F3. J) 

R I TE ( 6 , 70  ) 

70  FORMAT!////) 

C  PLOT  THE  IMPt'LS:  RESPONSE, 

CALLPLOT(U.O,-?.56,2) 

CALL°LOT( 0. 0, 0. 0, 2) 

00301=1, 128 

30  CALL3L0T(A< I ) / 4  J , -1*0,02, 2) 

CALL°LOT(0.0,-A,0,-3) 

C  3L0T  THE  MAGNITJDE  OF  TRANSFER  FUNCTION, 

00851=1, 128 
85  b< I ) =CAbS  <  H(  I  )  ) 


ooou  o  u  c>  o 
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x  M  s  o  .  0 
■J09UI  =  l,12fl 
90  Xin  =  MAXlF(BU>,X'1) 

CALL3LOT<0.0,-2.56,2> 

CALLpLOT(Q.O,O.U,2> 

00951=1,128 
b(  I  )=d(  1  )*2,5/X'i 
95  CALLPLOTtDC I >,-1*0.02,2) 

CALL°L0T(3. 0,4, 0,-3) 

iuo  continue 

CALLpLOT(-12,0,U,Q,-3) 

CALL°L0T(Q. 0,0. 0,999) 

CALLpLOTS(0,0,l) 

C  SHULATION  OF  PAYBACK. 

1>0200K  =  1,4 

mioi=i,i28 

110  1»(I)*(0. 0,0,0) 

sample  the  composite  array  to  retrieve 
all  the  samples  belonging  to  a  transfer  function, 
plot  the  magnitudes  of  the  sampled 
transfer  function, 

0012ni=r\,128,4 
120  0(I)=T(I) 

!>  01301  =  1, 128 
130  r ( I ) =C  AbS ( G (  I  )  ) 

X  H  a  U  ,  0 

i’0140I=l,  128 
140  Yii  =  MAXlF(H(  I  )  ,  XM) 

CALL°LOT<0.0,-2.56,2) 

CALL3LOT(0.0,0. J,2) 

CALLPLOT(0.0,-K»nt02,2) 

»0150I*K,12B,4 
"(  I  )=B( I )*2, 5/Xl 
OALLPLOT(B( I), -1*0, 02, 2) 

OALLpLOT(ti<  I),-<I*1)*0.02,2) 

OALL°LOT(0.0,-( 1*1)*0,02,2) 

150  CALLPLOT(0.0,-{ l+4)*0.02,2> 

CALLPLOT(0, 0,-4, 0,-3) 

FOJRlER  TRANSFORM  THE  SAMPLED  TRANSFER  FUNCTION. 
CALL  FOURIER  (G) 

3LOT  THE  OUTPUT  REPRESENT I NG  THE  PLAYBACK 
OF  THE  SYSTEM  WHEN  THE  TRANSFER  FUNCTIONS  ARE 
INDIVIDUALLY  ACCESSED, 

1-01601  =  1,128 
160  M  I ) =C  AdS ( G (  I  )  ) 

X  M  =  0 , 0 

U0170 I =1, 128 
170  XM=MAX1F(B( I >,X4> 
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1)01751*1,128 
175  »(!)*&(  I >*99/XM 

■RITE<6,  6 0 ) ( B ( 1 ) , I = 1 # 128 ) 

•'RITE (6#  70) 

CA|_L°L0T<  U.  0,-2, 56, 2) 

CALLPLOT(0.0,0. J,2) 

001801*1,128 

180  CALLPL0T(D( I)/4J,-I*0;02,2) 

CALL3L0T<3.0,4.U,-3> 

200  CONTINUE 

CALL»LOT(-12.0, J, 0,-3) 

CALLPL0T(G,C,99V) 

CALL=LOTS<0,0,1) 

plot  the  hagmitjde  of  the  composite  transfer  function 

ARRAY  REPRESENTING  THE  MULTIPLEXED  HOLOGRAM, 

D0210 1=1,128 
210  :*<  I  )=CAbS<T(  I  )  > 

XhrO.O 

f>0220I*l,l?8 
220  X  n  s  H  A  X 1 F  <  B  C  I  ),X*i) 

f'ALLPLCT{0.0,-2.56,2> 

CALL3LOT(&.0,O.U,2) 

1-02251=1,128 
■H  I  )=8(  l)*2,5/XM 
225  CALL3LOT(B(l)#-l*0,02.2) 

-ALLPLOT( 0 , 0, -4, 0, -3) 

FOJHIER  TRANSFORM  THE  COMPOSITE  ARRAY, 

CALL  FOURIER  <T> 

PLOT  THE  OUTPUT  REPRESENTING  THE  PLAYBACK 
OF  THE  SYSTEM  WHEN  ALL  THE  TRANSFER  FUNCTIONS 
“RE  SIMULTANEOUSLY  ACCESSED, 

CQ230 1=1, 128 
230  >'(  I  )=CASS(T<  I  )  ) 

X  K  =  0  ,  U 

>'02401  =1,128 
240  XHsHAXIF < B( I > , X«U 
002301=1,128 
250  !-  {  I  )sb(  I  )*R9/XM 

• R I T  E ( 6 ,  60)(B(i),I  =  l,128) 

CALL3LOT( 0 , 0, -2, 56, 2) 

CALL3LOT( 0 , 0, U . U , 2) 

1-0280  I  =1,128 

280  CALLPLOT(B< I )/4j,-J*0.02,2) 

CALL°LOT(0,0,4. J,-3) 

CALL  PLOTIO, 0,999) 

end 
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SUBROUTINE  FOURIER  (8) 

this  subroutine  calculates  the  fourier  transform 

OF  AN  ARRAY  OF  12B  ELEMENTS. 

THIS  PROGRAM  ALSO  SHIFTS  THE  ELEMENTS  IN  THE 
ARRAY  TO  TAKE  CARE  OF  THE  FFT  ALGORITHM  WHICH 
ASSUMES  THE  FIRST  ELEMENT  AS  THE  ORIGIN. 
U1MENSI0N  M ( 3 ) , S ( 32 ) * INV ( 32 )  ,8(128) 

TYPE  COMPLEXX.R 
OATA(M=7,0, 0) 

ITsO 

10  uO  20  1=1,64 
X  =  3<  I  ) 

4(1)=  8(1*64) 

M I*64)=X 
20  CONTINUE 

IF< IT. EO. 1)30,22 
\i  call  fft  aluoriihm  harm, 

22  -ALL  harm  (B,M, 1NV,S,1, IFfcRR) 

1  T  =  1 

Go  TO  10 
30  RETURN 
END 


APPENDIX  F 

Computer  Program  to  Generate  2-D  Composite 
Hologram  and  the  Output  of  a  Processor  Using  the 
Sampled  Transfer  Function  Approach 
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°  R  0  G  R  A  M  SAMPLER 

this  program  generates  and  plots  the  multiplexed 

TRANSFER  FUNCTION!  HOLOGRAM  OF  FOUR  IMPULSE 
RESPONSES  USING  THE  SAMPLING  TECHNIUUE  IN  THE 

frequency  plane,  the  program  also  simulates  and 

PLOTS  THE  SIMULTANEOUS  PLAYBACK  OF  ALL  THE 
TRANSFER  FUNCTIONS. 

DIMENSION  A(64,54>  ,3(32)  ,C(64) 

Tyofi  COMPLEX  A,X,CMPLX  ,B 
-RITE  (6 ,5>) 

5  K ORM AT (////) 

PU3. 1415926535 
DOlOOKsl, 2 
DOlOOLsl, 2 

READ  THE  IMPULSE  RESPONSE. 
REAU(5,10)((A(I,J),J=1,64>,I=1,64) 

10  F0RMAT(32(C(Fl, u ,  F 1 ,  0)  )  ) 

'•ULTIPLY  THE  IMPULSE  RESPONSE  WITH  a  CHECKERBOARD 
MASK  OF  *1  AND  -1  VALUES,  THIS  IS  PONE  TO  SPREAD  OUT 
The  TRANSFER  FUNCTION  MORE  EVENLY  IN  THE  FOURIER  PLaNE, 
(HOWEVER  THIS  WILL  NOT  AFFECT  THE  OUTPUT  OBSERVED 
AS  MAGNITUDE.) 

17  Doio  1=1,16 
001611*1,4 
D01B  0=1,16 
noiaojsi, 4 

18  A<<(I-1)*4*IJ),((J-1)*4+JJ)>=A(((I-1)*4+II), 

1< ( J-1)*4*JJJ )*(•!)**( I ♦ J ) 

OEVERATE  THE  TRANSFER  FUNCTION, 

Qallfqur I ER ( A  ) 

CONVERT  THE  VALUES  of  THE  TRANSFER  FUNCTION  INTO 
MAGNITUDE  AND  ANGLE,.. 

DQ50 I s 1, 64 
HOpUJsI, 64 
AM  =  CABS( A( I, J)  ) 

AA=CANG( A( I, J) ) 

IF(  AA)  20, 30, 30 
20  Aa=AA+2*PI 
30  AL=AM 

A( I, J)=CMPLX(AL, AA) 

50  CONTINUE 

sample  the  transfer  function  by  selecting  every 

ALTERNATE  ELEMENTS  IN  EACH  DIMENSION, 

STORE  THE  SAMPLES  ON  A  TaPE, 

DOSuIsK, 64, 2 

80  k'-RITE(3)(A(I,J),J  =  L,64,2) 

100  CONTINUE 
REWIND  3 


o  o  o  o  ci  o  o  o  rs  o  o 
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U  E  N  E  R  A  T  E  SAMPLED  COMPOSITE  TRANSFER  FUNCTION  HOLOGRAM. 
iJO200Ksl,  2 
302U  OLsl , 2 
IiO?0(UsK,64,2 

*EAD<3>y 
L'O2U0Jsl,32 
a ( I* ( ( J-l)*2+L) )=B( J) 

2U  0  CONTINUE 
HEJIMU3 

The  FOLLOWING  DJ  LOOP  GENERATES  THE  MULTIPLEXED 
hologram  plot  WHEN  the  INDEX  IS  1  and  SIMULATES 
the  simultaneous  playback  of  all  the  impulse 

RESPONSES  WHEN  INDEX  IS  2, 

210  Lq700Ls1*2 

scale  the  magnitudes  of  the  TRANSFER  FUNCTION 
to  a  value  OF  13, 

X  H  s  0 . 0 
302501=1,64 
i<02!?0  J  =  l,  o4 

250  XMsi14X1F<RFAL(A(  I,  J)),XM) 
uoson i si, 64 

C03un Jrl, 64 

AMs(REAL<A(I,J))/XM)*15,0 
3L0  All,  J>=CMPLX( AM, A  I  MAG ( A l 1 , J) >  > 

CALLPLOTS (0,0,1) 

PLOT  THE  ORIENTATION  MARKaR, 

:alL3LOT<  -G.5,'J.U,2  > 

CALLPLOTt  0 , 0  ,  0  ,  0 , 2  ) 

CALL°LOT (  0.0, -0.5, 2  ) 

CALLPLOT(  0.0, 0.0, 2  ) 

CaLLPLOT ( -1 , 5, 0 , 0 , -3) 

•■0*001*1,64 
•J0500Jsl,64 

IF  THE  PHASE  OF  THE  TRANSFER  FUNCTION  ELEMENT  IS 
BETWEEN  0  AND  120  DEGREES  RESOLVE  THE  VALUE  INTO 
C  COMPONENTS  ALONG  0  AND  120  DEGREES, 

Also 
*2  =  0 
A3  =  0 

a M  =  R E A L ( A ( I, J)  ) 

AA  =  AIMAG(A< I, J)  ) 

IF( AA.LT, 2*PI/3)330, 340 
330  Al*AM*(COS( A A )  +  (SI N( A A ) /SORT (3,0  )  ) ) 

A  A  s  2  *  P  I/3-AA 

A2=AM*(C0S<AA)*(SIN(AA)/SCRT(3,0))) 

G0TO370 


u  u  o  ou  o  u  u  o  u  u  o  ouo 
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IF  THE  PHASE  OF  THE  TRANSFER  FUNCTION  ELEMENT  IS 
HETHEEN  120  AND  240  DEGREES  RESOLVE  THE  VALUE 
INTO  COMPONENTS  along  120  AND  240  degrees. 

340  IF< AA.LT.4*PI/3)350,360 
350  AA=AA-2*PI/3 

A2=AM*(C0S(AA )*(SIN( AA)/S3RT(3. 0 ) > ) 

AAs2*PI/3-AA 

A3  =  AM*(C0S(  A  A  )  +  {SI'J(  AA)/SQRT(3,0  )  )  ) 

GUTO370 

IF  THE  PHASE  OF  THE  TRANSFER  FUNCTION  ELEMENT  IS 
BETWEEN  240  AND  360  DEGRFES  RESOLVE  THE  VALUE 
INTO  COMPONENTS  ALONG  240  AND  0  DEGREES, 

360  AA=AA-4*PI/3 

A3  =  AM*(COS(AA>4-{SIN(AA)/S3RT(3.0> ) ) 

AA  =  2*P  I/3-AA 

Al=AM*(COS(AA)+{SIN(AA)/SQRT(3,0))) 

••UANTIZE  THE  MAGNITUDE  OF  TRANSFER  FUNCTION  INTO  15 
STEPS  AND  PLOT  I  HE  RESOLVED  COMPONENTS  TO  A  WIDTH  OF 
0,05  INCHES  AND  THE  HEIGHT  PROPORTIONAL  TO  THE 
AGn  I  TULE  . 

370  IFIAl.LT.l, 0)40 j, 380 
360  N  M  s  A 1 

L'O390K  =  l  /  NM 

C ALLPLOT (-K*0, 01,-0.05,2) 

390  CALLPL0T(-K*0, 01, 0,0,2) 

400  CALLpLOT(U.O,-O.05,-3> 

IF (A2.LT ,l,0)43j,410 
410  ''  K  =  A  2 

U  0  4  2  0  K  *  1 ,  N  M 

0 A LLPL0T(-K*0 .01,-0.05,2) 

420  CALLPLOT(-K*0 .01,0,0,2) 

430  CALLPLOTIU .0,-0. 05, -3) 

IF (A3.LT .1.0)460,440 
440  J|v,sA3 

UOASOrsrl,  NM 

0 ALLPLOT (-K*0. 01,-0 .05, 2) 

450  CALL°LOT<-K*P.0l,0,0,2) 

f. 0 V £  THE  PEN  TO  THE  LOCATION  OF  THE  NEXT  ELEMENT, 

460  CALL=L0T( 0 , 0, -0 . 05, -3) 

500  CONTINUE 

.-‘OVE  THE  PEN  TO  THE  BEGlNlNG  OF  THE  NEXT  LINE, 
'ALLpLOT(-0,15,y,60,-3> 

600  CONTINUE 

RETURN  THE  PEN  TO  THE  STARTING  POSITION  AND  REMARK 

the  orientation  marker,  <  this  hill  check  the  total 

CUMULATIVE  ERROR  IN  THE  PLOTTER  POSITIONAL  ACCURACY , > 
CALLPLOTdl, 10,0, 0,-3) 


ci  rj  o  o  o 
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CaLLP|_0T<  0,0,-J,5»2  > 
callplok  o.o,o.n,2  ) 

CALL3L3T(  -Q.5,U,0,2  ) 

CALLPLOK  n.  0,0. 0/2  ) 

IF<L. £3.2)695, 6s3 

C  CONVERT  THE  MAGNITUDE  AND  PHASE  DF  THE  TRANSFER 

c  function  back  Into  real  and  Imaginary  parts, 

653  006541=1,64 
u0654J=l,64 
Am  =  REAL( A( I , J)  ) 

Aa*AImAG( A< I , J)  ) 

AX=AH*C0SF ( AA ) 

AYsAH*SINF( AA) 

654  A ( I, J)=CMPLX( AX, AY) 

FOURIER  TRANSFORM  THE  TRANSFER  FUNCTION  TO  SIMULATE 
THE  OUTPUT  OF  THE  SYSTEM  *hen  ALL  the  FUNCTIONS  ARE 
SIMULTANEOUSLY  jl ayed  sack, 

655  GALL  FOURIEK(A) 

CONVERT  ThE  OUTPUT  INTO  MAGNITUDE  AND  PHASE  AND 
RETURN  TO  THE  PLOT  ROUTINE  TO  PLOT  THE  OUTPUT, 
iJ  063  0  1=1,64 
D O 6 tt 0  J= 1 , 64 
AM  =  CAfaS( A( I, J)  ) 

AA=CAnG(A(I,J)> 

IF  (AA)660,6d0,6d0 
660  A  A  =  A A*  2*  P  I 
6ti0  A(  I,  J)=CMPLX(AM, AA) 

695  CALLPLOT(C), 0,999) 

700  CONTINUE 
=  ND 
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SURKOUTINE  FOUR  I FR  (A) 

r 

C  This  SUBROUTINE  COMPUTES  The  FOURIER  transform 
C  OF  A  2-U  ARRAY  3F  64  X  64  ELEMENTS. 

C  The  ohoGRaM  ALSu  SHIFTS  THE  QUADRANTS  OF  THE  ARRAY 

C  TO  TAKE  CARE  OF  THE  FFT  ALGORITHM  WHICH  ASSUMES  THE 

C  FIRST  ELEMENT  AS  THE  ORIGIN, 

C 

DIMENSION  a ( 64 , 64 ) , I  MV ( 16 ) , s ( 16 ) , M ( 3  > 

TYPE  COMPLEX  A,  x 
UATA(Ms6,6,0) 

I  T  =  0 
iFSETsl 
10  00201=1,32 
'>020  J  =  l,  32 
X=A( I , J) 

'(I,  J)=A<U32,J*32) 

A( 1*32, J*32)=X 
x  =  4< I , J*32) 

A ( I,  J«-32)=A< 1*32, J> 

20  A(J*32,J)=X 

IFUT.EQ.  1)40,30 

c  call  the  fast  fjurier  transform  algorithm  harm, 

3  0  CALLNARMf  A,M,  IN^S,  IFSET,  1F6R  n 
••H  I  TS  (  6#  35)  IFERR 
35  FuRhAT(5X,6HIFERPs, 13////) 

I  T  =  1 
L’OTGIO 
40  RETURN 
tND 
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APPENDIX  G 

Fabrication  of  2-D  Phase  Masks 


» 


» 


A  computer  program  to  plot  a  2-D  amplitude  mask  of  127  x 


127  elements  using  a  127  bit  Gold  code  given  in  Table  (2-2) 
has  been  described  [20].  Using  a  modified  version  of  this 
program  a  set  of  nine  plots  was  prepared  on  a  thick  bright 
white  drawing  sheet  to  a  size  of  3.81"  x  3.81"  each.  A 
typical  enlarged  plot  is  shown  in  the  Fig.  (G— 1 ) .  These 
nine  plots  were  mounted  on  a  white  poster  board  to  form  a 
3x3  array  with  center  to  center  distance  between  each  plot 
in  either  dimension  equal  to  13  inches.  This  array  of  plots 
was  then  photo-reduced  such  that  the  center  to  center  dis¬ 
tance  between  the  plots  is  equal  to  0.3  inches.  This  dimens¬ 
ion  was  chosen  to  match  a  fly's  eye  lens  array  with  which 
the  masks  were  to  be  used.  After  this  reduction  the  size 
of  each  cell  in  the  array  is  approximately  equal  to  18 
microns.  Thus  to  retain  good  resolution  after  photo  reduc¬ 
tion,  high  resolution  film  plates  type  Kodak  649F  was  used. 

The  exposure  details  and  the  processing  times  were  as  fol¬ 
lows  : 

Distance  between  the  camera  (fitted  with  a  50  mm  lens) 
and  the  plots:  94  inches. 

Exposure:  4  secs,  (plots  illuminated  by  diffused 

daylight  in  the  shadow  of  a  building  during  bright  sunlight) . 

Develop  in  Dll  solution:  12  minutes. 

Rinse  in  Kodak  stop  bath:  30  seconds. 

Rapid  fix  with  hardener:  5  minutes. 

Wash  in  running  water:  20  minutes. 

_  A 
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A  technique  to  fabricate  phase  masks  from  this  binary 
amplitude  mask  using  photo  resist  solution  has  been  des¬ 
cribed  [21].  It  was  observed  that  in  this  method  non  uni¬ 
form  thickness  of  the  coating  of  the  photo  resist  resulted 
in  variations  in  the  phase  difference  in  different  regions 
of  the  plate.  An  alternative  method  is  to  use  the  photo 
sensitive  property  of  gelatin  sensitized  by  a  dichromate 
solution.  A  number  of  papers  have  been  published  describ¬ 
ing  the  technique  of  fabricating  phase  holograms  using 
dichromated  gelatin  [22-26].  One  such  method  makes  use 
of  Kodak  649F  plates  in  the  starting  step  to  prepare  the 
sensitized  plates  [27].  The  advantage  of  using  the  649F 
plates  is  that  the  glass  base  of  the  film  plate  is  already 
coated  with  a  uniform  layer  of  gelatin  and  hence  the  prob¬ 
lem  of  coating  the  glass  with  a  uniform  gelatin  layer  is 
avoided.  The  detailed  processing  procedure  is  given  below: 
I.  Preparation  of  plates  coated  with  gelatin. 

(1)  Fix  a  64 9F  plate  in  rapid  fixer  with  hardener 
for  15  minutes. 

(2)  Wash  in  running  water  for  10  minutes. 

(3)  Soak  in  methyl  alcohol  for  10  minutes  with 
agitation. 

(4)  Soak  in  clean  methyl  alcohol  for  10  minutes 
with  agitation. 

(5)  Dry  in  a  vertical  position. 


» 
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I  At  the  end  of  these  steps  we  have  a  clear  glass  plate 

coated  with  a  layer  of  gelatin  on  one  side. 

II.  Sensitization  of  the  plate. 

(1)  Dissolve  10  gins  of  purified  ammonium  dich¬ 
romate  (Baker  brand  or  equivalent)  in  200  cc  of  dis¬ 
tilled  water.  Add  0.5  cc  of  photoflow  solution. 

(2)  Filter  the  solution. 

(3)  Place  the  glass  plate  with  the  gelatin 
side  up  in  a  flat  tray  and  pour  the  ammonium  dich¬ 
romate  solution  till  the  plate  is  completely  covered. 
Leave  it  in  this  position  for  5  minutes. 

(4)  Remove  from  the  solution  and  place  at  a 
small  inclination  (approximately  10°)  for  3  minutes 
to  let  the  excess  solution  to  flow  down.  Clean  the 
edge  of  the  plate  with  a  paper  towel. 

(5)  Place  in  a  light  tight  box  at  the  same  in¬ 
clination  as  in  step  4  above  for  24  hours. 

The  steps  3,  4  and  5  have  to  be  carried  out  under 
safelight  illumination  using  red  filter.  At  the  end  of 
these  steps  we  have  a  sensitized  and  prehardened  plate. 

III.  Exposure. 

(1)  Place  the  sensitized  plate  and  the  amplitude 
mask  such  that  their  emulsion  sides  are  facing  each 
other. 

(2)  Expose  for  12  minutes  under  a  500  watts 
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tungsten  filament  photo  lamp  placed  at  13  inches  dis¬ 
tance.  (Actually  several  exposures  ranging  from  4 
minutes  to  20  minutes  are  necessary  to  obtain  at 
least  one  mask  with  the  desired  phase  difference.) 

IV.  Development. 

(1)  Wash  in  clean  running  water  at  68°F  for  10 
minutes  under  safelight  illumination  using  red  filter. 

(2)  Soak  with  agitation  for  2  minutes  in  a  mix¬ 
ture  of  50%  isopropyl  alcohol  and  50%  water. 

(3)  Soak  with  agitation  for  10  minutes  in  100% 
isopropyl  alcohol. 

(4)  Pull  the  plates  out  of  the  alcohol  at  a 
rate  of  1  cm/min. ,  simultaneously  blowing  hot  air 
directed  at  the  surface  for  rapid  drying. 

These  steps  complete  the  process  and  a  phase  mask  is 
obtained. 

The  phase  mask  thus  fabricated  was  checked  in  a  Mach- 
Zehnder  interferometer  to  check  the  phase  difference  be¬ 
tween  the  elements.  However  in  view  of  the  extremely  small 
size  of  the  cells  it  is  difficult  to  resolve  the  fringe 
patterns  intersecting  each  cell.  Thus  a  reference  mark 
of  large  size  was  made  on  the  plots  and  was  used  as  ref¬ 
erence  to  check  the  phase  difference  between  the  exposed 
and  the  unexposed  parts.  However  because  of  the  following 
processing  problems,  phase  masks  to  the  desired  accuracy 
could  not  be  fabricated. 


(1)  Due  to  the  high  resolution  required  during  the 
fabrication  of  amplitude  masks  Kodak  649F  plates  were  used. 
But  these  high  resolution  plates  are  not  of  high  contrast 
type.  This  resulted  in  non  uniform  contrast  in  the  ampli¬ 
tude  mask  due  to  differences  in  illumination  and  hence 

non  uniform  exposure  of  the  dichromated  gelatin  plates. 

(2)  The  sensitivity  of  the  dichromated  gelatin  plate 
is  a  function  of  the  time  lag  after  sensitization.  Thus 
exposure  time  required  was  different  for  each  trial  depend¬ 
ing  on  this  prehardening. 

(3)  The  development  process  is  highly  sensitive  to 
the  temperature  and  the  Ph  of  water  [24]  and  it  was  not 
possible  to  control  these  parameters  in  the  existing  set 

Up. 

(4)  The  measurement  was  based  on  a  large  reference 
mark  made  on  the  plot.  But  due  to  difference  in  the  con¬ 
trast  of  this  reference  patch  with  those  of  the  actual 
cells  it  could  not  be  confirmed  whether  all  the  cells  in 
the  mask  have  the  same  phase  difference. 

All  these  problems  require  further  study  before  a 
standard  process  to  obtain  repeatable  results  may  be  fina¬ 


lized. 


